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We present transfer matrices for the zero-temperature partition function 
of the g-state Potts antiferromagnet (equivalently, the chromatic polynomial) 
on cyclic and Mobius strips of the square, triangular, and honeycomb lattices 
of width Ly and arbitrarily great length L^. We relate these results to our 
earlier exact solutions for square-lattice strips with Ly = 3, 4, 5, triangular- 
lattice strips with Ly = 2, 3, 4, and honeycomb-lattice strips with Ly = 2,3 and 
periodic or twisted periodic boundary conditions. We give a general expression 
for the chromatic polynomial of a Mobius strip of a lattice A and exact results 
for a subset of honeycomb-lattice transfer matrices, both of which are valid for 
arbitrary strip width Ly. New results are presented for the Ly = 5 strip of the 
triangular lattice and the Ly = A and Ly = 5 strips of the honeycomb lattice. 
Using these results and taking the infinite- length limit L^; — > oo, we determine 
the continuous accumulation locus of the zeros of the above partition function 
in the complex q plane, including the maximal real point of nonanalyticity of 
the degeneracy per site, W as a function of q. 
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I. INTRODUCTION 



The g-statc Potts antiferromagnet (AF) [1,2] exhibits nonzero ground state entropy, 5*0 > 
(without frustration) for sufficiently large g on a given lattice A or, more generally, on a 
graph G. This is equivalent to a ground state degeneracy per site W > 1, since 5*0 = ks In W . 
There is a close connection with graph theory here, since the zero-temperature partition 
function of the above-mentioned g-state Potts antiferromagnet on a graph G — G{V,E) 
defined by vertex and edge sets V and E satisfies 



where P{G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 



where n = |V"| denotes the number of vertices of the graph G and the symbol {G} formally 
denotes the set lim^i^oo G. The minimum number of colors that is necessary to color a graph 
G subject to this constraint is called the chromatic number of G, x{G). 

We represent a given strip as extending longitudinally in the x direction and transversely 
in the y direction, with width Ly vertices. Each strip involves a longitudinal repetition of m 
copies of a particular subgraph. For the square-lattice strips, this is a column of squares. It is 
convenient to represent the strip of the triangular lattice as obtained from the corresponding 
strip of the square lattice with additional diagonal edges connecting, say, the upper-left to 
lower-right vertices in each square. In both these cases, the length is L-^ = m vertices. We 
represent the strip of the honeycomb lattice in the form of bricks oriented horizontally. In 
this case, since there are two vertices in 1-1 correspondence with each horizontal side of a 
brick, = 2m vertices. 

The general structure of the Potts model partition function for cyclic strips of the square 
lattice of width Ly, as a sum of powers of eigenvalues of a formal transfer matrix multiplied 
by certain coefficients c^'^\ < d < Ly, was given in Refs. [7] (see also [8]). The present 
authors (unaware of this finding in [7]) rediscovered the result in Ref. [9] and showed that it 
applies also to cyclic strips of the triangular and honeycomb lattices [10]. The coefficients 
are polynomials of degree d in the variable q given by [7,9] 



Z{G,q,T^O)pAF^P{G,q) 



(1.1) 



[3]- [6]. Thus 



W{{G},q)^ hm P{G,q) 



l/n 



(1.2) 




j=0 



(1.3) 
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with Un{x) being the Chebyshev polynomial of the second kind. The first few of these 
coefficients are c(°^ = 1, c^^^ = g — 1, c^^^ = — 3g + 1, and c^^^ = — bq^ + <oq — l. The c^'^^'s 
play a role analogous to multiplicities of eigenvalues XA,Ly,d,j, although this identification is 
formal, since c*^"'-' may be zero or negative for the physical values q = 1, 2, 3. For example, as 
shown in eqs. (2.19), (2.20) of Ref. [9], if g = 2, then c^"^) = -1 for = 2 or 3 mod 4, and if 
g = 3, then c^"^ = — 1 for = 3 or 5 mod 6 and c'-'^^ = —2 for ci = 4 mod 6. In general, c^'^^ 
vanishes at g = qd,k, where qd,k = 4 cos^(7r/c/(2(i + 1)), k = 1,2, so that c^'^^ is positive 
for g > 4 for arbitrary d. 

This structure also holds for the chromatic polynomial, with the difference that, for 
Ly >2, the number of eigenvalues for each d is smaller than the number for the full partition 
function. We determined this number for strips of the square and triangular lattice in Ref. 
[9], for strips of the honeycomb lattice in Ref. [10], and for self-dual strips of the square 
lattice in Ref. [11] (where we also noted how these numbers fill out the entries in the relevant 
BrateUi diagrams) . The chromatic polynomial for a cyclic strip of the regular lattice A has 
the form 

Ly np{A,Ly,d) 

P(A,L,xm,cyc.,g)=5^cW J] (AA,L„.,,r (1-4) 

d=0 j=l 

Let G' — (y, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E. Then P(G, g) can be written as the sum [12,13] 

P(G,g)= ^g'^(^')(-l)l^'l (1.5) 

G'CG 

where k{G') denotes the number of connected components of G' and \E'\ denotes the number 
of edges in the set E'. Since we only consider connected graphs G, we have k{G) = 1. The 
formula (1.5) enables one to generalize g from Z+ to M+. The zeros of P{G,q) in the 
complex g plane are called chromatic zeros. We denote the continuous accumulation set of 
these zeros in the n — oo limit as B, which is the continuous locus of points where M^({G}, g) 
is nonanalytic. The maximal value of g where B intersects the (positive) real axis is labelled 
1c{{G}). This locus occurs as a solution to the degeneracy in magnitude of A's of maximal 
magnitude [14]- [17]. 

In this paper we present transfer matrices Ta,Lj^ for strips of the square (sq), triangu- 
lar (tri), and honeycomb (he) lattices of fixed transverse widths Ly and arbitrarily great 
length with periodic (cyclic) longitudinal boundary conditions. We also give results for the 
corresponding strips with twisted periodic (Mobius) longitudinal boundary conditions. We 
relate these results to our earlier exact solutions for square-lattice strips with Ly — 3,4,5; 
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triangular-lattice strips with Ly = 2,3,4, and honeycomb-lattice strips with Ly = 2,3 and 
periodic or twisted periodic boundary conditions. We give a general expression for the chro- 
matic polynomial of a Mobius strip of a lattice A which is valid for arbitrary strip width Ly. 
New results are presented for the Ly = 5 strip of the triangular lattice and the Ly = 4 and 
Ly = 5 strips of the honeycomb lattice. Using these results and taking the infinite-length 
limit Lx oo, we determine certain properties of the continuous accumulation locus of the 
zeros of the above partition function in the complex q plane, including the maximal point 
of nonanalyticity of the degeneracy per site, 1^ as a function of q. We also give a general 
trace formula for the honeycomb strip which is valid for arbitrary Ly. Chromatic numbers 
for these lattice strips are given in the appendix. 

There are several motivations for this work. One is that, as noted above, the Potts 
antiferromagnet has the interesting property of nonzero ground state entropy for sufficiently 
large q, which is an exception to the third law of thermodynamics [18,19]. Via eq. (1.2), exact 
calculations of q) yield calculations of I^IG}, q) and thus give insight into this property. 
The constraint that no two adjacent vertices have the same value of q becomes less and less 
important as q ^ oo, and in this limit, q) — > As q decreases, q) remains 

a real analytic function of q down to the point qc{{G}), where it is nonanalytic. In earlier 
work, we studied how qc{{G}) depends on the boundary conditions used for various lattice 
strips [20]- [39]. It was found that a convenient property of strips with periodic (or twisted 
periodic) longitudinal boundary conditions was that they always defined a value of qc{{G}), 
and, furthermore, for a given type of lattice, this point was observed to be a monotonically 
nondecreasing function of the strip width for all cases calculated. By carrying out exact 
calculations of P{G,q) and determining the analytic structure of W{{G},q), one can study 
how qc{{G}) varies with increasing strip width and how it approaches the values expected for 
the corresponding two-dimensional lattices, namely qdsq) = 3 [40] and qdtri) = 4 [41,42]. 
This interpolation property is an especially interesting use of these exact results since the 
resultant values of qc{{G}) interpolate between the value qc — 2 for one-dimension and 
values on two-dimensional lattices. Similar results were found for the ground state entropy 
itself [20,44-47]. This should be contrasted with other properties; for example, the critical 
temperature of the Potts ferromagnet on an infinite-length lattice strip is zero for any width 
Ly regardless of how great. 

For the honeycomb lattice, formal arguments give qdhc) = (3 -|- V5)/2 ~ 2.618 [48,43]. 
The honeycomb value is only formal because if q is not a positive integer, then one cannot 
use the expression of the partition function for the Potts antiferromagnet in terms of a sum 
of positive-definite Boltzmann factors and instead must use the formula [12] Z{G, q, v) — 
S^G'CG?'^^'^''*'^'^''' where v — — 1 with K — 13 J, where (3 — {ksT)'^ and J is the spin- 
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spin coupling (which reduces to eq. (1.5) for the T = antiferromagnet, where K ~ — oo). 
However, for —1 < v < and non-integral g, this does not, in general, define a Gibbs measure 
[48,49]. This non-integral property of qdhc) for the two-dimensional honeycomb lattice and 
the associated subtleties make it especially useful to have explicit calculations of qc values 
for finite-width strips of this lattice, to check that these arc consistent with the expected 
behavior, according to which they approach the above value of qdhc) as Ly gets large. 

A second motivation related to this is that when one generalizes q from a positive integer 
to a complex variable, one sees that the point qc{{G}) is the maximal point where a certain 
nonanalytic boundary B crosses the real axis (discussed further below); thus, another purpose 
of these exact calculations is to get further insight about this locus B. As is well known, 
slight changes in coefficients of polynomials can have drastic changes on the positions of the 
zeros of the polynomials, which means that one must have the exact expression for P{G, q) 
to study these zeros and their accumulation set B. Since it was found that for (the L^, — > oo 
limit of) strips with free longitudinal boundary conditions, B does not necessarily cross the 
real axis, so that no qc{{G} , q) is defined [20,23], one cannot carry out this study in the same 
manner for these strips (although for sufficiently wide strips with free longitudinal boundary 
conditions, arcs on B have endpoints that are often close enough to the real axis to allow 
one to extrapolate and define an effective qc). 

A third motivation is that certain structural features of the transfer matrices can allow 
one to obtain general formulas that are applicable for arbitrarily large strip widths. We 
have already used this feature in Ref. [39], where we presented transfer matrices whose 
eigenvalues are \\^Ly,d with degree d = Ly — 1 loi several lattice strips with periodic boundary 
conditions, including those of the square and triangular lattices. Thus another purpose of 
exhibiting explicit transfer matrices for < d < Ly — 2 here is to make them available for 
a wider community, leading, hopefully, to some advances beyond those wc have made with 
d = Ly, Ly — 1 in the construction of explicit formulas for \K,Ly,d valid for arbitrary Ly. All of 
these degree-d sectors contribute to the partition function (here, the chromatic polynomial), 
so that they are of physical interest, although an elementary result is that the degeneracy 
per site W only depends on the (maximal eigenvalue in the) d — Q sector. 

The chromatic polynomial for the cyclic and Mobius strips of the square lattice were 
calculated (i) for L^^ = 2 in [50] (see also [51]); (n) for Ly — ?, cyclic in [29,30], and Mobius 
case in [31]; (ni) for Ly — A cyclic/Mobius in [35], (iv) for Ly — 5 cyclic/Mobius in [39]. 
The chromatic polynomials for the triangular lattice strips were calculated for (i) Ly — 2, 
cycfic/ Mobius in [30], (ii) = 3, 4 in [36]; and for the honeycomb lattice with (i) Ly — 2, 
cychc/Mobius in [27], (ii) Lj, = 3 in [10]. (Refs. [27,28] showed that the structure (1.4) holds 
more generally for cyclic strip graphs composed of iterated subgraphs that are not parts of 
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regular lattices; we shall not need this degree of generality here.) We have proved that for 
these strips the set of eigenvalues are the same for cyclic (eye.) and Mobius (Mb.) boundary 
conditions and have given the transformation rules for how the coefficients change for strips 
of the square lattice when one changes from cyclic to Mobius boundary conditions [9]. 

In our previous works on chromatic polynomials we have sometimes given transfer ma- 
trices (e.g. [39]) but often have expressed our results for the \A^Ly,d,j in terms of generating 
functions or solutions to algebraic equations. In one respect this is a maximally compact 
way of presenting the results, since for a given set of np{A, Ly,d) AA,Lj,,d,i's, it suffices to 
give the eoeffieients of the algebraic equation of degree np{A, Ly,d) and for the generating 
function, it suffices to give 2np{A, Ly,d) coefficients, which are polynomials in q. Since the 
generating function is a rational function in gr, it has the appeal that one always deals with 
polynomials. For the expression of the chromatic polynomial in terms of transfer matrices, 
since Ta^l^^^ has dimension np(A, Ly, d) x np{A, Ly, d), one must specify a greater number of 
polynomials, namely np{A, Ly,dy polynomials. As Ly increases, the number of expressions 
that one must give thus grows quadratically rather than linearly, as in the more compact ap- 
proach using generating functions or algebraic equations for the \K,Ly,d- However, offsetting 
this disadvantage, the transfer matrix method has the advantage that some of the entries are 
rather simple, and, moreover, one can sometimes spot useful patterns in the matrices that 
expedite or confirm the construction of general formulas that are applicable for arbitrarily 
large strip widths Ly. 

We mention some previous work. Matrix methods and related recursive linear algebraic 
techniques for calculating chromatic polynomials were used in early work [50-52] and more 
recent papers [53]- [55]. The coloring matrix technique of Ref. [52] was applied in Refs. 
[44]- [46] to obtain upper and lower bounds on W for various (infinite) two-dimensional 
lattices. This application relied upon the fact that these coloring matrices are non-negative, 
so that one could use the Perron- Frobenius theorem in analyzing the eigenvalues. Related 
recursive linear algebraic methods were used in Refs. [23-25] and in our subsequent papers 
on chromatic polynomials. Transfer matrices for chromatic polynomials were used in Refs. 
[41,42,56]. In Ref. [57] transfer matrices for the Potts model were developed and were apphed 
there and in Refs. [58,59] to calculate chromatic polynomials for strips of the square and 
triangular lattices with free longitudinal boundary conditions and free or periodic transverse 
boundary conditions. These have been termed transfer matrices in the Fortuin-Kasteleyn 
representation (see eq. (1.5)). In Refs. [60,61], transfer matrix methods were used to calculate 
full Potts model partition functions with arbitrary q and v for strips of the square and 
triangular lattices with free longitudinal boundary conditions. Using transfer matrix methods 
together with the sieve methods of Refs. [53]- [55] (see also [62-64]), we have proved theorems 
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that determined T\^Ly,d ior d = Ly — 1 and arbitrarily large Ly for several lattice strips with 
periodic boundary conditions [39]. 

Some remarks are in order here concerning an important respect in which the transfer 
matrices presented here differ from the conventional transfer matrices in statistical mechanics 
and coloring matrices in mathematical graph theory. Consider a strip graph, and let the 
spins on each transverse slice of this strip be labeled a^^y where the longitudinal position x 
is fixed and the transverse position y varies. As conventionally defined, the transfer matrix 
of the (zero-field) Potts model at some temperature T is defined as the matrix T^^x+i = 
{Wx,y}\e~'^'^\{o-x+i,y}) , where -J Y^i^ij) ^aiaj, c^i = 1, Q are the spin variables on each 
vertex and (ij) denotes pairs of adjacent vertices. Since there are q^^ spin configurations on 
each shce, this is a q^'« x q^y dimensional matrix. Although this dimension depends on q, 
the individual elements themselves do not explicitly depend on q. For physical temperatures 
< T < oo this matrix is non-negative, and at nonzero temperatures it is positive-definite. 
For T > 0, on a strip of finite width, for which this is a finite-dimensional matrix, one can then 
make powerful use of the Perron-Frobenius theorem on positive-definite matrices to conclude 
that there is a unique real positive eigenvalue of maximal magnitude. Indeed, this theorem 
provides a standard way to prove that the free energy of the Potts model or other spin model 
with short-range interactions on a strip of infinite length and finite width is analytic at all 
finite temperatures. Similarly, the coloring matrices used in Refs. [52,44-46] for studies of 
chromatic polynomials and the degeneracy per site W on various lattices are non-negative 
matrices, and this property was necessary for the upper and lower bounds derived in these 
papers. In contrast, the transfer matrices in the Fortuin-Kasteleyn representation are not 
positive-definite, and hence the Perron-Frobenius does not apply to them. The fact that 
it does not apply is crucial for the property that the locus B crosses the positive real axis 
and for the existence of qc{{G}), since this locus occurs where there is a switch between 
at least two distinct eigenvalues of maximal magnitude and hence non-uniqueness of such 
eigenvalues. 

From the exact expressions for P(G, g), one can evaluate eq. (1.2) to get W and analyze 
the accumulation locus B of zeros of P{G, q) in the complex q plane. We have done this 
for cyclic/ Mobius strips in a number of previous works. These loci are different for strips 
with different boundary conditions, e.g., free, cylindrical, cyclic/ Mobius, and toroidal. As 
we have shown, since the A's are the same for the cyclic and Mobius strips of a given lattice, 
it follows that the locus B is the same for the infinite-length limits of cyclic and Mobius 
strips of a given lattice. In our previous work we have found a number of general properties 
of the locus B, including the property that for strips of regular lattices with periodic or 
twisted periodic longitudinal boundary conditions the locus B is comprised of closed curves 
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that enclose various regions and pass through g = and at a maximal real point qc which 
depends on the lattice, as well as possible other intermediate points. The point qc{{G}) is 
important since it separates the interval q > qc{{G}) on the positive real q axis where the 
Potts model (with q extended from Z+ to M) exhibits nonzero ground state entropy (which 
increases with q, asymptotically approaching 5*0 = feeing for large q) from the interval 
< g < qc{{G}) in which 5*0 has a different analytic form. 



II. GENERAL STRUCTURE 



Our results for chromatic polynomials are obtained as special cases, for the zero- 
temperature antiferromagnet, of general transfer matrices that we have calculated for the full 
temperature-dependent Potts model partition function on cyclic and Mobius lattice strips 
[65]. The latter results, and the general method, will be reported in a companion paper [66]. 
The motivations for presenting the special cases of these results for the zero-temperature 
Potts antiferromagnet (chromatic polynomial) have been given in the introduction. In ad- 
dition to these physics motivations, a relevant point is that the transfer matrices for the 
chromatic polynomials are subtantially smaller in dimension and simpler in structure, de- 
pending only one one variable instead of two, than the transfer matrices for the full Potts 
model partition function. Another reason for separating these analyses of the full Potts 
model and the chromatic polynomial has to do with the determinants det{Tz,K,Ly,d) and 
det{Tp^\^Ly^d)', in both the case of the full Potts model and the chromatic polynomial, the 
individual eigenvalues Xx,A.Ly,d for most strip widths are not expressible in explicit alge- 
braic form because the corresponding characteristic polynomials are of fifth order or higher. 
However, we have found a most interesting property, that the products of the eigenvalues 
that define the determinant det{Tz,K,Ly,d) have a very simple form. When one specializes to 
V — —1, ioT all values of Ly and d = ...Ly — 1, except the lowest case, Ly = 1, some of 
these eigenvalues vanish, so that det{Tz^A^Ly,d) = at v = —1, reflecting the fact that the 
transfer matrices Tp^j^^Ly,d are of smaller dimension than Tz,K,Ly,d- But in this v — —1 case, 
when one removes the zero columns and corresponding rows of Tz,K,Ly,d to form the trans- 
fer matrix Tp^a^l^^^ for the chromatic polynomial, the resultant determinant is, in general, 
nonzero. And again we find a very interesting feature, namely that although the eigenvalues 
Xp,A,Ly,d themselves are roots of characteristic polynomials that are often of fifth or higher 
order, precluding solutions in terms of explicit algebraic expressions, the products that define 
dei(Tp^A,Lj,,(i) are often quite simple, especially for the strips of the triangular and honeycomb 
lattices. This is another motivation for our presenting results for the chromatic polynomials 
separately from those for the full Potts model partition function. 
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The chromatic polynomial for cyclic strips can be written in the form [7,9] 

Ly 

P{A, Ly X m, eye, q) = Y, c^'^TrHT^^Ly^T) (2.1) 

where Tf^^i^^d is the transfer matrix, with eigenvalues \A^Ly,d- We denote the characteristic 
polynomial of T\^Ly,d in the variable z as CP(T/y^Ly,d, z). We have shown that the full transfer 
matrix has a block structure formally specified by 

Ly 

TA,Ly = ®\[TA,Ly,d (2.2) 

where the product n^A.-t-j/.d means a set of square blocks, each of dimension, c^'^\ of the 
form XA,Ly,d,j times the identity matrix. The dimension of the total transfer matrix, i.e., the 
total number of eigenvalues XA,Ly,d,j in eq. (1.4), counting multiplicities, is thus 

Ly Ly 

dim{TA,Ly) = dMTA,Ly,d) = Ly, d) . (2.3) 

d=0 d=0 

Each matrix TA,Ly,d is defined relative to a basis of d-color assignments to a path graph 
with Ly vertices corresponding to a transverse slice across the strip. Here one introduces 
partitions of the vertices in these path graphs. Using these as bases, we calculate the transfer 
matrix for each level (= degree) d, as discussed in our Rcf. [66]. For the chromatic polyno- 
mials of the square and triangular lattices, we have the additional requirement that adjacent 
vertices cannot connect to each other. We list graphically all the possible partitions for 
Ly = 2,3,4 strips in Figs. 1 to 3, where white circles are the original Ly vertices and each 
black circle corresponds to a specific color assignment. We denote these partitions VLy,d for 
2 < Ly < 5 as follows: 

7^2,0 = {1} , 7^2,1 = {2; 1} , 7^2,2 = {1, 2} (2.4) 

P3,o = {/;13}, P3,i = {3;2;l;T3} 

7^3,2 = {2, 3; 1, 3; 1, 2} , Vs,s = {1, 2, 3} (2.5) 

P4,o = {/; 13; 14; 24} , ^4,1 = {4; 3; 2; 1; 13, 4; l3; Ti; 24; 24, 1} 
P4,2 = {3, 4; 2, 4; 1, 4; 2, 3; I, 3; 1, 2; T3, 4; 1, 24} 
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P4,3 = {2, 3, 4; 1,3, 4; 1,2, 4; 1,2, 3} , 



P4,4 = {i,2,3,4} 



(2.6) 



^5,0 = {/; 13; 14; 15; 24; 25; 35; 24, 15; 135} 

= {5; 4; 3; 2; 1; 13, 5; 13, 4; T3; 14, 5; l4; T5; 24, 5; 24; 24, 1; 25; 25, 1; 35; 35, 2; 35, 1; 
24,T5;I35} 

^5,2 = {4, 5; 3, 5; 2, 5; 1, 5; 3, 4; 2, 4; 1, 4; 2, 3; 1, 3; 1, 2; 13, 4, 5; T3, 5; l3, 4; l4, 5; 24, 5; 24, 1, 5; 

1,24; 1,25; 2, 35; 1,35; 35,1,2} 
7^5,3 = {3, 4, 5; 2, 4, 5; 1, 4, 5; 2, 3, 5; 1, 3, 5; 1, 2, 5; 2, 3, 4; 1, 3, 4; 1, 2, 4; 1, 2, 3; T3, 4, 5; 1, 24, 5; 

1,2,35} 

n,4 = {2, 3, 4, 5; 1, 3, 4, 5; 1, 2, 4, 5; 1, 2, 3, 5; 1, 2, 3, 4} 

n,5 = {1,2,3,4,5} (2.7) 



where partitions are separated by a colon, and each overhne corresponds to a color assign- 
ment. That is, for d = there is no explicit color assignments; for d = 1, one vertex, or 
one connected set of vertices has its color specified; for d = 2, two vertices or two separate 
connected sets of vertices have their colors specified; and so forth for higher values of d. 
Therefore, in eqs. (2.1) to (2.4), there is no overline for d — 0, one overline for each partition 
oi d— 1, two over lines for each partition oi d — 2, etc. 

The number of partitions np{A, Ly, d) is the dimension of the transfer matrix T\^L^d, 



In Ref. [9] we determined the number np{A, Ly,d) (labelled simply as np{Ly,d)) for the 
lattices A = sq,tri). Some special cases for A = sq,tri are [9] as follows; analogous results 
for the honeycomb lattice were given in Ref. [10]. 



dim{TA,Ly,d) = np{A, Ly, d) . 



(2.8) 



np{A,Ly,Ly) = 1 



(2.9) 



np{A, Ly, Ly -1) = Ly 



(2.10) 
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np(A, Ly, 0) = np{A, Ly- 1,1)^ Ml,-i 
where M„ is the Motzkin number in combinatorics, given by 



(2.11) 



(2.12) 



where C„ 



1)^^(^^) is the Catalan number. (We use the same symbol C„ for the 
circuit graph; the meaning will be clear from context). For reference, the first few Motzkin 
numbers are M„ = 1, 2, 4, 9, 21, 51 for n = 1, 6, and we use the formal definition Mq = 1. 

Note that the number np{tri, Ly,0) is the same as the dimension of the corresponding 
transfer matrix for the triangular-lattice strip with free boundary conditions obtained in 
Ref. [57], while the transfer matrix for the free strip of the square lattice has a dimension 
(for which a general formula was given in Ref. [60]) which is smaller than np{sq, Ly,0) for 
Ly > A. Consistent with this, for Ly > 4, we find that the characteristic polynomials for 
Tsq^o and Tgq^^fi factorize into parts such that one factor is the characteristic polynomial for 
the transfer matrix of the corresponding free strip, as will be seen below. 
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FIG. 2. Partitions for the Ly = 3 strip. 
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FIG. 3. Partitions for the Ly = A strip. 



As in earlier work, we define Np^A,Ly,\ as the total number of distinct eigenvalues of T\^Ly, 
i.e. the sum of the dimensions of the submatrices Ta l d, modulo the multiplicity c^^^ : 



(2.13) 



d=0 



Since we will use these results here, we include the relevant values in the Appendix. For 
cyclic or Mobius strips of the square and triangular lattices. 



{Ly-J) 



3=0 



{mLy-2j)\ 



for A = sq, tri 



(2.14) 



where [z/] denotes the integral part of u. 

For the corresponding strips of the honeycomb lattice, we calculated np{hc, Ly,d) and 
Np,hc,Ly,x in Ref. [10]. The results that are needed here are hsted in the appendix. These 
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correspond to the numbers of partitions which allow the connections between vertices 2 and 
3, between vertices 4 and 5, etc. Therefore, in addition to the partitions given in eq. (2.5), 
the other partitions for the Ly — 3 strip of the honeycomb lattice are 

Vhc,3,o = {23} , 7^^0,3,1 = {23; 23, 1} , Vhc,3,2 = {1, 23} . (2.15) 

In addition to the partitions given in eq. (2.6), the other partitions for the Ly = 4 strip of 
the honeycomb lattice are 

VhcA,o = {23; 23, 14} , Vhc,4,i = {23, 4; 23; 23, 1; 23, 11} 

Vhc,4,2 = {23, 4; 23, 1, 4; 1, 23} , VhcA,3 = {1, 23, 4} 

In addition to the partitions given in eq. (2.7), the other partitions for the Ly 
the honeycomb lattice are 

Phc,5,o = {23; 45; 13, 45; 23, 14; 23, 15; 23, 45; 145; 235; 245; 23, 145} 
Phc,5,i = {23, 5; 23, 4; 23; 23, 1; 45; 45, 3; 45, 2; 45, 1; 13, 45; 45, 13; 14, 23, 5; 23, 14; 23, T5; 

23, 45; 45, 23; 23, 45, 1; 145; 235; 235, 1; 245; 245, 1; 23, 145} 
P/^c,5,2 = {23, 4, 5; 23, 5; 23, 1, 5; 23, 4; 23, 1, 4; I, 23; 3, 45; 2, 45; 1, 45; 45, 2, 3; 45, 1, 3; 45, 1, 2; 

13, 45; 23, 14, 5; 23, 45; 23, 1, 45; 45, 1, 23; 1, 235; 1, 245} 

Vhc,5,3 = {23, 4, 5; 23, i, 4, 5; 1, 23, 5; 1, 23, 4; 2, 3, 45; 1, 3, 45; 1, 2, 45; 45, i, 2, 3; 1, 23, 45} 

7'/.c,5,4 = {l,23,4,5;l,2,3,45} . (2.17) 
For cyclic strips of the square and triangular lattices, general properties include first [9] 

Ly 

dim{TA^Ly) = P{^, Ly x m, eye, q)m=o = c^'^'^np{A, Ly, d) 

d=0 

= PiTreeLy, q) = qiq - 1)^^"' , A^sq, tri (2.18) 
where TrecL^^ denotes a tree graph (here a path graph) with Ly vertices; and 

Ly 

Tr{TA,Ly) = P{A, Ly X m, eye, q)m=i = J2 c^^^Tr{TA,Ly,d) = , A = sq,tri. (2.19) 

d=0 

12 



(2.16) 
5 strip of 



For cyclic strips of the honeycomb lattice [10] 



Ly 

dim{Thc,Ly) = P{hc, Ly x m, eye, q)m=o = ^ c^'^^np{hc, Ly, d) 

d=0 

for Ly odd 

^ ^ (2.20) 

— 1)^ ^ for Lj^ even 

Here wc give another general result for these strips of the honeycomb lattice, governing the 
trace of the transfer matrix: 

Ly 

Tr{Thc,Ly) = P{hc,Ly x m,cyc.,q)m=i = ^c^'^^Tr{Thc,Ly,d) 

d=0 

= P(Tree2L„ q) = q{q - l)'^-"^ . (2.21) 

This is proved by the same coloring matrix methods that we used in Ref. [9] and [10]. The 
reason for the difference relative to the corresponding formula (2.19) for the cychc strips of 
the square and triangular lattices is that for m = 1 on those strips the longitudinal edges 
connect each vertex to itself, so the chromatic polynomial vanishes, but, in contrast, for 
the strip of the honeycomb lattice, the longitudinal edges have interior vertices of degree 
2 associated with multiple edges. Because of this, and taking into account the elementary 
theorem that the chromatic polynomial of a graph remains unchanged if one replaces any 
edge by multiple copies of this edge, it follows that the resulting coloring is described by 
the chromatic polynomial P(Tree2Li;, Q')- For the reader's convenience in checking that our 
results satisfy these general trace formulas, we shall list traces of our transfer matrices below. 
In general, 

Ly 

dem,Ly) = Y[[dem,Ly,d)f' ■ (2.22) 

We next discuss Mobius strips. One of the major results of the present paper is a 
general formula, eq. (2.26) together with eqs. (2.23)-(2.25), for the partition function of the 
zero-temperature g-state Potts antiferromagnet on the Mobius strip of a regular lattice A of 
arbitrary width as well as arbitrary length. This formula involves a set of Ta ^ which occur 
once in each m-fold product (with the Ta ^ thus raised to the power m — 1) and which 
differ from the Ta.Lj^,^ in the interchange of certain columns, encoding the Mobius property 
of reversed-orientation periodic longitudinal boundary conditions. In our previous papers 
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involving Mobius strips of the square and triangular lattices, we expressed the chromatic 
polynomial in a form analogous to (1.4) but with different sets of coefficients. In Ref. [9] we 
derived a set of transformations that governed how the coefficients for a certain d change 
when one changes the longitudinal boundary condition from cyclic to Mobius for the square- 
lattice strip. However, as noted there, this did not apply to the triangular-lattice strip; for 
that strip, if one expresses the chromatic polynomial in the same manner for the cyclic and 
Mobius cases, as a sum of m'th powers of eigenvalues Xtri,Ly,d,j, then some of the coefficients 
become algebraic, rather than polynomial, functions of q. This was evident in the simplest 
nontrivial example, the Ly = 2 strip [30]. The general formalism that we introduce here for 
Mobius strips has the advantage that the coefficients stay in the set of c^'^^'s for triangular 
(and honeycomb) as well as square lattice strips. For Ly > 3, the TA^Ly,d matrices have 
eigenvalues that, in general, are not just related to those of T\^Ly^d by possible sign changes. 
In the limit of infinite length, m — > oo, because only the T\^Ly,d matrices are raised to a 
power that goes to infinity, it is only the eigenvalues of these matrices that determine locus 
B. Thus, the result of our previous papers, that the locus B is the same for the infinite-length 
limit of the strip of the lattice A for cychc or Mobius boundary conditions, is again evident 
in our present formulation. 

We proceed to the details. For Mobius strips one set of horizontal edges reverses the 
order of connection. This corresponds to the exchange of the pair of bases which switch 
to each other when the vertices reverse the order, i.e., the set of bases which do not have 
self-reficction symmetry with respect to the center of the tree with Ly vertices. For example, 
among the partitions for the Ly = 2 strip in Fig. 1, the first partition 2 and the second 
partition 1 in 7^2,1 must be exchanged under this reflection. The pairs of partitions for the 
Ly = 3 strip in Fig. 2 are the flrst partition 3 and the third partition 1 in V^^i, and the flrst 
partition 2, 3 and the third partition 1, 2 in P3.2. This corresponds to the exchange of these 
pairs of columns of T^^L^^d, and these matrices will be denoted as Ta^Lj,,^ [66]. For general 
Ly we have the following changes of coefficients for the square, triangular and honeycomb 
lattices [9]: 



(2.23) 



1 < A; < 



- L 2 J 



(2.24) 



,(2fe+l) (fc+l) 



< A; < 



1 



(2.25) 



L 2 
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We thus obtain the important general formula, in terms of transfer matrices, for the chro- 
matic polynomial for a Mobius strip of the lattice A, which is valid for arbitrary width 
Ly-. 

P{KLy X m,Mh.,q) = c^'^TT[{T^,L^,or~'f^,Lj 

[{Ly-l)/2] 
[Ly/2\ 

- c('^-'^rr[(rA,L„2.)"^-^Tk,i^,2.] . (2.26) 

d=l 

(where the factor c*^*^^ = 1 is displayed explicitly for uniformity). 

For the square lattice, the eigenvalues of Tsq^Ly,d are the same as those of Tsq^Ly,d except 
for possible changes of signs. The number of eigenvalues with sign changes is the number 
of column-exchanges from Tsq,Ly4 to Tsq^Ly,d- Denote the number of eigenvalues that arc the 
same for Tsq ^^^^ and Tsq^Ly,d as np{sq, Ly, d, +), and the number of eigenvalues with different 
signs as np{sq, Ly, d, —). It is clear that 

np{sq, Ly, d, +) + np{sq, Ly, d, -) = np{sq, Ly, d) . (2.27) 

Define 

Anp{sq, Ly, d) = np{sq, Ly, d, +) - np{sq, Ly, d, -) (2.28) 

which gives the number of partitions having self-reflection symmetry. For example, among 
the partitions for the Ly — 2 strip in Fig. 1, the partition I in 7^2,0 and the partition 1,2 
in 7^2,2 have self-reflection symmetry. Among the partitions for the Ly — 3 strip in Fig. 2, 
it includes the two partitions in Vsfi, the second partition 2 and the fourth partition 13 in 
7^3,1, the second partition 1,3 in 7'3_2, and the partition in Vs^s- We list Anp{sq, Ly,d) for 
1 < Ly < 10 in Tabic III. The total number of these partitions for each Ly will be denoted as 
ANp^Ly, and we have AA^p2n-i = A.Np^2n for n > 0. The relations between Anp{sq, Ly,d) 
are 

Anp{sq, 2n — 1,0) = Anp{sq, 2n, 0) , for < n 
Anp{sq, 2n — 1, 2m — 1) = Anp{sq, 2n, 2m) , for 1 < m < n 
Anp{sq, 2n — 1, 2m) — Anp{sq, 2n, 2m — 1) , for 1 < m < n 



15 



Anp{sq, 2n +1,0) = 2Anp{sq, 2n, 0) + Anp{sq, 2n, 1) , for < n 



Anp{sq, 2n + l,m) = Anp{sq, 2n, m — 1) + Anp{sq, 2n, m) + Anp{sq, 2n, m + 1) , 

for < m < 2n + 1 
Anp{sq, 2n + 1, 2m) = Anp{sq, 2n + 1, 2m + 1) , for < m < n 
Anp{sq, 2n, 0) = Anp{sq, 2n, 2) 

Anp{sq, 2n, 2m - 1) = Anp{sq, 2n, 2m + 2) , for < m < n - 1 . (2.29) 

Wc also list np{sq,Ly,d,-\~) and np{sq,Ly,d,—) for 2 < Lj, < 10 in Table IV. Notice that 
np{sq, Ly, 0, +) is the number of \p^sq,FF,Ly proved in [60]. Recall the numbers of Xp^sq,Ly,j for 
the Mobius strips of the square lattice with coefficients ic*^^', defined as np^Mb{Ly, d,±) = 
np,Mb{sq, Ly,d,±), have been given in [9]. With the eqs. (2.23) to (2.25), the relations 
between np{sq, Ly, d, ±) and np^Mb{sq, Ly, d, ±) are 

np,Mb{sq, Ly, 0, ±) = np{sq, Ly, 0, ±) + np{sq, Ly, 2, +) 

np^Mb{sq,Ly,k,±) = np{sq,Ly,2k - 1,±) +np{sq,Ly,2k + 2,^) , l<k< ^ ^ (2.30) 

We illustrate the application of this general formalism in the section where we give results for 
specific transfer matrices. For Mobius strips of the square lattice, the sum of the coefficients 
is given by 

[iLy-l)/2] 

P{sq,LyXm,Mb.,q)m=o = c^°^Anp{sq,Ly,0)+ c^'^^^^ Anp{sq, Ly,2d + 1) 

d=0 

[Ly/2] 

- J2 c^'^~^^Anp{sq,Ly,2d) 

d=l 

= c(°) [Anp{sq, Ly, 0) - Anp{sq, Ly, 2)] 

KLy + l)/2] 

+ J2 ^^'^^ [^np{sq, Ly, 2d -I)- Anp{sq, Ly, 2d + 2)] 

d=l 

for Ly even 

(2.31) 

P(Tree(L^+i)/2, g) for Ly odd 
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which agrees with Theorem 5 in [9] because of eqn. (2.30). This m = equation also holds 
for the triangular lattice strips. If wc take m = 1 in eqn. (2.26), the chromatic polynomial 
is zero for the square lattice with odd Ly or for the triangular lattice with any Ly. This is 
what we expect since these Mobius strips have at least one edge connect a vertex to itself. 

For the honeycomb lattice with Ly even, the magnitudes of eigenvalues of Thc,Ly,d and 
Thc,Ly,d ai'c again the same. We denote the number of eigenvalues that arc the same as 
np{hc, Ly, d, +), and the number of eigenvalues with different signs as np{hc, Ly, d, —). Sim- 
ilarly to eqns. (2.27) and (2.28), we have 

np{hc, Ly, d) — np{hc, Ly, d, +) + np{hc, Ly, d, — ) 

Anp(/ic, Ly, d) = np{hc, Ly, d, +) — np{hc, Ly, d, — ) . (2.32) 

The values of Anp{hc, Ly,d) for even Ly up to Ly = 10 arc listed in Table V. The values 
Anp{hc, Ly, d) can be obtained from Anp{sq, Ly, d). For example. 



Anp{hc,2,d) 


= Anp{sq, 2, d) 


Anp{hc,4:,d) 


— Anp{sq, 3, d) + Anp{sq, 4, d) 


Anp{hc, 6, d) 


— Anp{sq, 4, d) + Anp{sq, 6, d) 


Anp{hc, 8,d) 


= Anp{sq, 5, d) + Anp{sq, 6, d) + Anp{sq, 7, d) + Anp{sq, 8, d) 


Anp{hc, 10,0?) 


= Anp{sq, 6, d) + 2Anp{sq, 8, d) + Anp{sq, 10, d) 



The relations between Anp{hc, Ly, d) are 

Anp{hc, An, 0) = AAnp{hc, An -2,0) + 2Anp{hc, An - 2, 1) , for < n 

Anp{hc, An, 2m — 1) = np{hc. An, 2m) 

= 2{Anp{hc, An - 2, 2m - 1) + Anp{hc, An - 2, 2m)) 

+Anp{hc, An - 2, 2m - 2) + Anp{hc, An - 2, 2m + 1) , 

for 1 < m < 2n 
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Anp{hc, 4n + 2, 0) = Anp{hc, 4n + 2, 2) 



= Anp{hc, 4n, 0) + Anp{hc, 4n, 2) , for 1 < n 
Anp{hc, 4n + 2, 2m - 1) = np{hc, 4n + 2, 2m + 2) 

= Anp{hc, in, 2m — 1) + Anp{hc, in, 2m + 2) , for 1 < m < 2n 

(2.34) 

np{hc, Ly, d, +) and np{hc, Ly, d, — ) for even 2 < Ly < 10 are listed in Table VI. For Mobius 
strips of the honeycomb lattice with even Ly, the sum of the coefficients is given by 

[iLy~l)/2] 

P{hc,LyXm,Mb.,q),ri=o = c''^^Anp{hc,Ly,0)+ ^ c^'^^^'^ Anp{hc, Ly,2d + 1) 

d=0 

[Ly/2] 

- c^'^~^^^rip{hc,Ly,2d) 

d=l 

= c(°)[Anp(/ic, Ly, 0) - Anp{hc, Ly, 2)] 

KLy + l)/2] 

+ J2 c^"*^ [Anp{hc, Ly,2d-1)- Anp{hc, Ly, 2d + 2)] 

d=l 

for Ly/2 odd 

/ ^Ly/4 (2-35) 

— 1)1 for I/j//2 even 

It will be convenient to introduce the following polynomial: 

Dn - = ^ ') (2-36) 

where P{Cn, q) — {q — 1)" + — 1)(— 1)" is the chromatic polynomial for the circuit graph 
with n vertices, C„. Various properties of these polynomials were given in Refs. [25,26]. The 
first few DnS are Di = 0, D2 = 1, = q — 2, D^ — q^ — 3gf + 3. In addition to the identities 
given in Refs. [25,26], another identity is 

Dn + = (g - 1)"-' ,n>2. (2.37) 

Certain linear combinations of these polynomials also occur in the entries of the transfer 
matrices. These can be seen to arise as powers of basic linear factors such as g — a plus 
integer constants; e.g., q^ — 4q-\-b — {q — 2Y -\- 1, q^ — Qq-\- — {q — 3)^ + 1, etc. 
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We shall write out the entries of the transfer matrices explicitly if this can be done in 
the space available. In cases where this is not feasible, we shall express the entries in terms 
of the above polynomials, together with the shorthand notation 

Qa^q-a, qa,b ^ (q - a)iq - b) (2.38) 

and additional shorthand defined in the last appendix, including the polynomials F^^n, G^^n, 
and others. 

III. SOME TRANSFER MATRIX RESULTS FOR ARBITRARY Ly 

A. d = Ly 

For arbitrary Ly, eq. (2.9) shows that for d — Ly the transfer matrix T/^,L^,d=L is one- 
dimensional, i.e., a scalar. We have previously proved that [39] 

XA,Ly,Ly = {-lf' for A = sq,tri. (3.1) 

Our results also led to the inference that [27,10] 

Aa,l„l, = 1 for A = hc. (3.2) 

These results can also be seen to follow immediately from the transfer matrix method. 

B. d = Ly-l 
1. Square lattice 

For arbitrary Ly, the dimension of the transfer matrix is specified by eq. (2.10) as 
np{A, Ly, Ly — 1) = Ly for A = sq,tri. For this case we have given the transfer matrices 
TK,Ly,d=Ly-i for cychc strips of the square lattice in eqs. (7.1.1) and (7.1.4)-(7.1.6) of Ref. 
[39] and for cyclic strips of the triangular lattice in eqs. (7.3.1)-(7.3.4) of Ref. [39]. General 
formulas for eigenvalues and determinants were also given in Ref. [39]. Since we will need 
these here, we list them below: 

K,Ly4=Ly~i,j = {-lf^^\q - as,,Ly,j) ,l<J<Ly (3.3) 

where 

as,,Ly,, = 1 + 4 cos^ ,l<J<Ly. (3.4) 
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Hence, the determinant and trace of Tsq,Ly,Ly-i are 

Ly 

det{Tsg,Ly,Ly-i) = ~ (^sg,Ly,j) (3.5) 

Tr{T,,,,y,Ly-i) = (-1)^^+M2 + Ly{q - 3)] . (3.6) 

Although our general result (3.3) exhibits an explicit factorization of the characteristic poly- 
nomial of Tsq^Ly,Ly-i into linear factors, this involves, in general, transcendental numbers 
(trigonometric functions) as coefficients. A few of these are of the form ±(g — k) where k is 
an integer; e.g., for Ly — 2, d — 1, one has Xsq^2,i,j — ^ ~ Q,^ ~ Q and for Ly — 3, d — 2, one 
has Xsq,3,2,j = ?~l)?~2,g — 4. Some others involve terms of the form ±(g — a) where a is 
an algebraic number; e.g., for Ly — A, d — 3, one has A4^3j = 1 — q,3 — q,3 ± y/2 — q. In 
Table VII we indicate the factorization properties of the various characteristic polynomials 
of Tsq^Ly,d in factors with integer coefficients; for the illustrative case Ly = 4, d = 3, this is 
CP{Tsg,4,3,z) = {z-l + q){z-3 + q){z^-2{3-q)z + q'^-6q + 7), symbolized as (1^,2) in the 
table. We observe that the factorizations often correspond to the numbers np{sq, Ly,d,+) 
and np{sq, Ly, d, — ) in Table IV, although sometimes there are additional factorizations. 

2. Triangular lattice 
For Ttri,Ly,Ly-i we proved that [39] 

det{Ttri,Ly,Ly-i) = (? - 2)^(g - 3)^^-^ (3.7) 

Tr{Ttr^,Ly,Ly-l) = (-1)^^+M3 + Ly{q - 4)] . (3.8) 

Since a matrix and its transpose have the same eigenvalues, it is a convention which one 
is presented. The transfer matrices obtained directly from our calculation of the analogous 
matrices for the full Potts model are the transposes of those given in [39]. This just amounts 
to a different convention for the ordering of the basis of configurational states. 

3. Honeycomb lattice 

Here we give an exact determination of Thc,Ly,d ior d — Ly — 1 which is valid for arbitrary 
Ly. This goes beyond our results in Ref. [39], which presented such formulas for the square 
and triangular lattices. We first observe that this matrix has dimension 
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for L„ odd 



dim{Thc,Ly,Ly-i) = np{hc, Ly, Ly - 1) = (3.9) 

^Ly — 1 for Ly even 

This can be derived in a straightforward manner from results in [10]. It is convenient to work 
with two auxihary sets of matrices, each of dimension nz{hc, Ly, Ly — 1) = 2Ly — 1, which 
are special cases of the transfer matrices for the full Potts model on strips of the honeycomb 
lattice, Shc,Ly,Ly-i,j for j — 1, 2. The elements of these matrices will be given below. We 
then form the product 

^hC,Ly,Ly — l — Shl.,Ly,Ly-l,2ShC,Ly,Ly-l,l ' (3.10) 

In this matrix certain columns are identically zero. The number of such zero columns is 

Ly-l 



for L„ odd 



2 ^ y 

nz{hc, Ly, Ly-l) - np{hc, Ly, Ly - I) = \ (3-11) 

^ for Ly even 

The transfer matrix Thc,Ly,Ly-i is then the submatrix of Shc,Ly,Ly-i formed by removing these 
columns with vanishing entries and the corresponding rows, which operation we denote with 
a subscript red. for "reduced" : 

Thc,Ly,Ly — \ — Sfic,Ly,Ly—l,red. ■ (3-12) 

To determine Thc,Ly,Ly-i, it thus suffices to give the Shc,Ly,Ly-i,j matrices for j — 1,2. For 
the relevant range Ly >2 where the strips of the honeycomb lattice are defined we calculate 

{Shc,Ly,Ly-i,i)j,j^2-q for l<j<Ly-l (3.13) 
(Q \ _/2-g for even 

[^hC,Ly,Ly-l,l)Ly,Ly " ^ 1 _ ^ fof odd ^"^-^^^ 

{Shc,Ly,Ly-l,l)2j-l,2j = {Shc,Ly,Ly-l,l)2j,2j-l = 1 for l<j< [Ly/2] (3.15) 

iShC,Ly,Ly-l,l)Ly+2j-l,Ly+2J^1^0 l<j<[Ly/2] (3.16) 

{SHc,Ly,Ly-l,l)Ly+2J,Ly+2J^l ^OT l<j<[{Ly-l)/2] (3.17) 

{Shc,Ly,Ly-l,l)2j-l,Ly+2j-l = hc,L y ,L y -I ^l) 2j ,L y+2j -1 ^OT 1 < j < [Ly/2] (3.18) 
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{ShC,Ly,Ly-l,l)2j,Ly+2j = {S hC, Ly , Ly - 1 ,l) 2j + 1 , L y + 2j = "1 1 <j < [{Ly " l) / 2] 

{Shc,Ly,Ly-l,l)Ly+2j-l,2j-l ^ {S hc,Ly ,Ly -1 ,l) Ly +2j -1 ,2j ^ "l ^OT 1 <j < [Ly/2\ 

{ShC,Ly,Ly -1,2) 1,1 = 1 — ? 

{Shc,Ly,Ly-i,2)j,j = 2-q for 2 < j < Lj, - 1 
, , ( 2 — q for Lj^ odd 

{^hC,Ly,Ly-l,2)Ly,Ly = ^ 1 _ ^ fof Lj, eVCIl 

iShc,Ly,Ly-l,2)2j,2j+l = ('S'fec,Lj„Lj,-l,l)2j+l,2j = 1 for 1 < J < [{Ly - l)/2] 
('S'^c,Lj„Ly-l,2)Ly+2j-l,Ly+2j-l = 1 for 1 < j < [Ly/2] 
{Sf,,,Ly,Ly-l,2)Ly+2j,Ly+2j^0 for 1 < j < [{Ly - 1) /2] 
{Shc,Ly,Ly-l,2)2j-l,Ly+2j-l = {Shc,Ly,Ly-l,2)2j,Ly+2j-l = -1 for 1 < i < [-^j//2] 
{Shc,Ly,Ly~l,2)2j,Ly+2j = ('S'/ic,Ly,Ly-l,2)2j+l,Ly+2j = for 1 < j < [{Ly - l)/2] 
{Shc,Ly,Ly-l,2)Ly+2j,2j = ('S'fec,Ly,Ly-l,l)Ly+2j,2j+l = "1 for 1 < J < [{Ly - l)/2] 

with all other elements equal to zero. Using our general formulas, we find that 
det{Thc,Ly,Ly-i) = {q- for odd Ly>3. 

For even Ly, det{Thc,Ly,Ly-i) is proportional to {q - 1)^^F4,3 = {q- l)^y{q^ -Aq + 5). 

We next proceed to list our explicit calculations of transfer matrices. Since we have 
completely determined the transfer matrices in the cases d = Ly, where it is a scalar, and 
for the case o? = L^^ — 1, we do not list the former for each specific Ly. 
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IV. STRIPS OF THE SQUARE LATTICE 

A. Ly = 2 

We recall first that for Ly — 1, an elementary calculation yields P{sq, 1 x m, eye, q) — 
{q - 1)"* + {q - 1)(-1)"* so Tsq^ifi = g - 1. The chromatic polynomials P{sq, 2 x m, BC, q) 
for cyclic and Mobius EC's were given in Ref. [50] One has np(A, 2, 0) = 1, np(A, 2, 1) = 2, 
np(A2, 2) = 1 for A = sq (or tri), Tsq^2,o = - 3g + 3, and 

T.,..=-{'':,\-_\) (4.1) 

yielding the eigenvalues \sq,2,i,i = 1 — 9 and Asg^2,i,i = 3 — g. The determinant det{Tsq^2,i) = 
(g — l)(g — 3) and trace Tr{Tsq^2,i) = — 2(g — 2) follow as the Ly = 2 special case of our 
results in [39], viz., eqs. (3.5) and (3.6). For the present strip, the general methods discussed 
above yield, for the Mobius strip, the matrix Tsq^2,i which has the first and second columns 
interchanged. 

B. Ly = 3 

The chromatic polynomials P(sg, 3 x m, BC, g) were given in [29,30] for cychc EC's and 
in [31] for Mobius EC's. We find 

/(g-2)(g2-3g + 4) g2-4g + 5\ . . 

Ts,,3,o - ^ g - 2 J ^^-^^ 



sq,3,l 



/-(g2-4g + 5) g-2 
g-2 -(g-2)2 
-1 g-2 

V 1 1 



-1 
g-2 



2-g\ 
1 

-(g2-4g + 5) 2-g 
1 g-2y 



(4.3) 




T.,,3,2= -1 g-3 -1 (4.4) 



Although Tsq^3fi is the same as the analogous transfer matrix for the free strip of width 
Ly = 3, Tgq^Lyfi is larger than the free-strip transfer matrix for Ly>A. Note that the upper 
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left-hand 3x3 block of T^q,3,i is symmetric. The matrix Tgq^s^i is obtained by exchanging 
columns 1 and 3 of 2^5^,3,1, and 2^5 3 2 is obtained by exchanging columns 1 and 3 of ^55 3 2. 
We have 

det{Ts,,3,o) = {q- - + 13q - 11) (4.5) 
det{T,g,s,i) = -{q - l){q - 2)\q' - 9q^ + 29q' - 4Qq + 22) (4.6) 



det{T,,,s,2) = (g - l)(g -2){q- 4) (4.7) 
rr(r,,,3,o) = (g-2)(g^-3g + 5) (4.8) 

Tr(T,,,3,i) = - 13g + 16) . (4.9) 

Tr{Tsq^3^2) is given by the Ly — 3,d — 2 special case of the general formula (7.1.28) of [39], 
given above as eq. (3.6). The characteristic polynomial CP{Tsq^2,,i, factorizes into polyno- 
mials of degree 1 and 3 in z. This is indicated in Table VII. The eigenvalue corresponding 
to the linear factor is Asq^3^i^i = —{q — 2)^ [29]. 



C. Ly = A 



The chromatic polynomials P{sq, 4 x m, BC, q) for BC=cyclic and Mobius were given in 
[35]. We find 



( S21 {q-2)pe rso (g - 2)p6 \ 

q-2 {q-2f 3-q 1 

-1 2-q g^-5g + 7 2-q 

\q-2 1 3-q {q-2fj 



where the polynomials S21, P&i and r2o are defined in the appendix. 
We have 

det{Tsq,i,o) = (g - l)(g - 3){q^ - IQq^ + 112g^ - 449g^ + 1130g^ - 1829g^ 



(4.10) 



+1858g^ - 1084g + 279) 



(4.11) 



and 
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Tr{T, 



sqAfij 



= - + 24g2 - 45g + 36 . 



(4.12) 



As indicated in Table VII, the characteristic polynomial of Tsq,4fl consists of factors of degree 
1 and 3 in z. The eigenvalue which is the root of the linear factor is \sq,i,o,i = (? ~ !)(? ~ 3); 
the eigenvalues of the cubic factor are the same as the eigenvalues for the free strip of width 
Ly — A determined in Ref . [23] . The matrix 3^9,4,0 is obtained by exchanging columns 2 and 
4 of Tsg,4,o- 



■59,4,1 



/ 


-ri3 


-^4,3 


-q2 


1 


—^4,3 


12 


—^4,3 


-ll 


12 \ 




-^4,3 


— ?2-^4,3 




-12 


12 


-ll 


53 


12 


-1 








— Q'2-^4,3 


-^4,3 


-1 


12 


Iz 


-ll 


12 




1 


-q2 


-^4,3 


-^^13 


12 


-ll 




12 


~G^4,3 




-1 











-12 
















-1 


Q2 


?2 


12 


-12 


ll 


-g3 


1 


1 




-1 


-1 


-1 


-1 


-12 


-12 


^4,3 


-12 


-12 




Q2 


Q2 


12 


-1 


1 


1 


-1z 


ll 


-12 


V 











-1 














-12 1 



(4.13) 



The shorthand notation used in eq. (4.13) and below, such as Fm,n and Gm,n, is defined in 
the last appendix. Note that the upper left-hand 4x4 block of Tgq^^^i is symmetric. We 
calculate 



c?et(T,,,4,i) = {l- lf{l - ^f{l^ - 3g + 3)2 
X {q^ - 17g^ + 125g^ - 520g^ + 1342g^ - 2206g^ + 2261^^ - 1325g + 341) 

Tr{T,q,4,,) = -4g3 + 27q^ - 69q + 65 . 



(4.14) 
(4.15) 



The matrix Tsg,4,i is obtained by exchanging columns 1 and 4, 2 and 3, 5 and 9, and 6 and 8 
of Tsq,4,i. It is ironic that although the chromatic polynomial is a special case of the full Potts 
model partition function, the determinants of the transfer matrices TA,Lj,,d for the degree-d 
subspaces for these square-lattice strips, such as eq. (4.14), are more complicated than the 
simple expression that we found for the determinant of the transfer matrix Tz,sq,Ly,d in Ref. 
[66], 



det{Tz, 



sq,Ly,d) 



V 



Ly [riz {Ly ,d)-nz(Ly-l,d)] 



nz(Ly-l,d) 



(4.16) 



where here nz{Ly,d) = nz{A, Ly,d) for A = sq,tri. Furthermore, the expression for 
det{Tsq^Ly,d) for these square-lattice strips gets more complicated as Ly increases, in contrast 
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to eq. (4.16) and also in contrast to the results that we find for det{TA,Ly,d) for A — tri, he 
(see below). 

As indicated in Table VII, the characteristic polynomial of Tsq^4^i consists of factors of 
degree 4 and 5. 



sq,4,2 



( ^4,3 


2 


- q 


1 










(1 


- 2 




\ 


2-q 




5g + 


6 3-q 


2 


-q 1 







-1 







1 


3 


- Q 


Ps 




1 3-g 


1 


Q 


- 2 


Q 


- 2 





2 


- Q 


1 




-2)2 2-q 




















1 


3-g 


2 


- Q Q2,3 


2-q 









-1 










1 




2-q 


-^4,3 







Q 


- 2 


-1 




-1 


-1 










2 


- Q 







V 







-1 




-1 


-1 







2 


-9/ 



(4.17) 



where Fm,n and pg are defined in the last appendix. The upper left-hand 6x6 block of T,^ 4 2 
is symmetric. The matrix T^g 4 2 is obtained by exchanging columns 1 and 6, 2 and 5, and 7 
and 8 of Tsq^4^2- CP{Tsq^4^2, z) has factors of degree 3 and 5. We have 



det{Tsq,4,2) = {q- lf{q - 3f{q^ - 7?' + 15? - 11) 
X (g^ - 16?^ + 106?^ - 378g^ + 788g^ - 967?^ + 653g - 189) 

Tr(T,,,4,2) = - 29q + 38 



(4.18) 

(4.19) 



- sg,4,3 



/g-2-1 0\ 

-1 g-3 ~1 
-1 q~3 -1 
\0 -\ q-2 ) 



(4.20) 



Prom the general formulas (3.5) and (3.6), we have ^6^(75^,4,3) = (g — l)(g — 3)(g2 — 6g-|-7) 
and rr(r,,,4,3) = -2(2g-5). 



D. = 5 



The chromatic polynomials P{sq, 5 x m, BC, q) for BC=cyclic and Mobius were given in 
[39]. In addition to the general determination of Tsq^Ly,Ly-i, the transfer matrix for Tsq^^^^ 
was given in [39]. We display here our calculation of Tsq^^fi. 
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/ ^4,3^^11 


S31 


S47 




9ip7 


S47 


S31 


^19 


92P8 


-^4,3 


?2-^4,3 


-?2,3 


2g-5 


92 


93 


92 


1 


9i 






92^4^3 


-P8 


-4 


93 


-1 


-92 


-92 


1 


(12 


~G'4^3 


^^20 


92 


— <54,3 


92 




9f 


<i 


?2 


-92,3 


2^3 


9i 


-92,3 


92 


92,3 


1 


-q2 


-1 


93 


-P8 


-9i 


92^^4,3 


-9i 


-92 


-92 


-^4,3 


?2 


93 


2q-5 


92 


-92,3 


92-^4,3 


1 


9i 











1 











92 





V -1 


-q2 


-93 


5-2g 


1 


-93 


-92 


-1 


-9i 



(4.21) 

where the various p, r, s polynomials are defined in the appendix. The characteristic poly- 
nomial of Tsqfifi consists of factors of degree 2 and 7 in z. The eigenvalues that are the roots 
of the degree-7 factor are the same as those found in Ref. [23] for the corresponding Ly = 5 
free strip of the square lattice. 

The transfer matrices T^^ for d = 1,2 are both of dimension 21 x 21 and are too lengthy 
to present here; they are available from the authors. We find 

det{Tsg,5,o) = -(9 - 1)^(9 - 2f{q^ - 13g^ + 69g^ - 191g^ + 292g^ - 236g + 79) 
x(q^^ - 34g^^ + 538g^^ - 5259?^^ + 35541?^^ - 176036?^° 
+660682g9 - 1914798g^ + 4324155g^ - 7615130g^ + 10381339g^ 

-10768339g^ + 8235159^2 - 4388527?^ + 1459163g - 228580) . (4.22) 

The determinant det{Tsq^5,i) factorizes as (g — 1)^(9 — 2)^P2o-P29, where P20 and P29 are 
polynomials in q of degree 20 and 29 in q. det{Tsq^5^2) factorizes as {q — 1)^(9 — 2)Pi9P23, 
where P19 and P23 are polynomials in q of degree 19 and 23. For d = 3 we have 

det{T,q,5,3) = (9 - 1)'(9 - 2) (9' - 20g^ + 171q^ - 817g^ + 2389g^ 
-4387g^ + 4954g2 - 3155g + 869) (g" - 29g^° + 375g^ - 2853g^ + 14188g^ 

-48439g^ + 115934g^ - 194762g^ + 225461g^ - 171683g^ + 77617g - 15835) (4.23) 
and, from eq. (3.5), 
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det{Ts,,,,4) = (g - l){q^ -5q + 5){q^ - 7q + 11) (4.24) 
Wc obtain the traces 

Tr{Tsq,5,o) ^q^ - + A2q^ - 120g^ + 195? - 137 (4.25) 

Tr{Tsq,r^,i) = -5q^ + 46g^ - 180g^ + 346g - 269 (4.26) 

Tr{Ts,,5,2) = W - 74g2 + 198g - 189 (4.27) 

rr(r,,,5,3) = -W + 51g - 69 (4.28) 
with Tr(Tgg^5^4) being given by eq. (3.6). 

V. STRIPS OF THE TRIANGULAR LATTICE 

A. Ly = 2 

The chromatic polynomials P{tri, 2 x m, BC, q) for BC=cyclic and Mobius were given 
in [30]. One has Xtri,2,o = {q ~ 2)^, and, as a special case of our general result in eqs. (7.1.1) 
and (7.3.1)-(7.3.4) of [39], 

r„,. = -f^-^!-n (6.1) 



-1 q-2 



with eigenvalues 



1 r 



5-2q±y/9-Aq (5.2) 



Note that we use a different ordering convention for the basis configurations here than the 
one wc used in Ref. [39], so that Ttri^Ly,Ly-i in terms of our present basis is the transpose 
of the corresponding matrix given in [39]. This has no effect on the eigenvalues, since the 
eigenvalues of a matrix A and its transpose are the same. We have det{Ttri^2,i) = (Q'~2)^ 
and Tr{Ttri,2,i) = 5 — 2q, in accordance with the general formulas (3.7) and (3.8). For the 
Mobius strip, the Ttri^2,i matrix is obtained from Tfri,2,i by interchanging the two columns, 
as was true for the square lattice strip. However, in contrast to the case of the square-lattice 
strip, where this just reverses the sign of one of the eigenvalues, here the Ttri,2,i matrix has 
different eigenvalues than Ttri,2,i (namely, (l/2)(3 — g± i/5g^ — 22q + 25 )). It is instructive 
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to display how our general formula operates. For comparison, we first display the explicit 
chromatic polynomial for the cyclic strip [30]: 



P{tri, 2 X m, eye, q) = c^^\Xtri,2,or + c^'^ (At„,2,i 



.(1) 



(A. 



tri,2,l,2 



+ c' 



(2) 



(5.3) 



If one expresses the chromatic polynomial for the Mobius strip as a sum of powers of the 
same eigenvalues as for the cyclic strip, the result is [30] 

(?-!)(? -3) 



P{trt,2 X m,Mb.,q) = c^^\Xtri,2,oy 



(AtH,2,l,l)"* - (Atri,2,l,2)' 



.(0) 



(5.4) 



One sees that the coefficient of the second and third terms is no longer in the set of c^^^'s 
and, indeed, is not a polynomial function of q. Of course, the square root in the denominator 
of this coefficient cancels so that, as must be true, the chromatic polynomial is a polynomial 
in q. As discussed in [31], this can be seen as a consequence of the theorem on symmetric 
functions of roots of a polynomial equation whose coefficients are polynomials in q; although 
these roots are not, in general polynomial functions of q, the symmetric functions can be 
expressed in terms of polynomials in q. In the present case, the second and third terms 
can be written in a manifestly symmetric manner by observing that the denominator of the 
coefficient is equal to A^^i, 2,1,1 — Atri,2,i,2- With our general formula (2.26) we express the 
chromatic polynomial in a form that keeps the coefficients in the set of c^^^ polynomials. 



B. Ly = 3 



The chromatic polynomials P{tri, 3 x m, BC, q) for BC=cyclic and Mobius were given in 
[36]. We have 



- tri,3,0 



(g-2)(g2_55 + 7) {q-?>f 
2-q q-3 



(5.5) 



This yields the same eigenvalues that were earlier calculated for the free Ly = 3 strip in 
Ref. [23] and can be taken to be the transfer matrix for that strip. 



T, 



tri,3,l — 



f -{q''-6q + 10) -(g-2)(g-3) q-2 3-q\ 

q-3 -{q-2){q- 3) -{q - 2)(g - 3) 3 - g 

-1 q-2 -(q-2){q-3) 3-q 

\ 2 2-q 2-q q-3 J 



(5.6) 



The characteristic polynomial CP{Ttri,3,Q, z) is a quadratic, while CP{Ttri,3,i, z) has a 
hnear and cubic factor, as indicated in Table VIII. We note that for odd Ly the factorizations 
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given in Table VIII correspond to the numbers np{sq, Ly, d, +) and np{sq, Ly, d, — ) given in 
Table IV, but there are no nontrivial factorizations for even Ly. 



As a special case of our general formulas given in eqs. (7.3.1)-(7.3.4) of [39], we have 

'-tri,3,2= I -J- q-^ q-z. \ (5.7) 

We calculate 



q-3 q-4 q-2 
-1 q-A q-2 
-1 q-2 



detiTtr^,3,o) = iq-2fiq-3) 



(5i 



det{Ttri,3,i)^-{q-2f{q-3f 



det{Ttri,s,2) = (g - 2f{q - 3) 



Tr{Ttri,s,o) =q'- Iq^ + 18g - 17 



Tr{T,ri,3,i) = -3g2 + I7q - 25 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



Tr{Ttri,3,2) = 3(g - 3) . 



(5.13) 



C. Ly = 4 



The chromatic polynomials P{tri, 4 x m, eye, q) were given in [36]. The transfer matrices 
that we give below, together with our general procedure for obtaining T^^^d from Tj^^l^^, yields 
a new solution for P{tri, 4 x m, Mh., q). We find 



Ttri,4,0 — 



({q-2){q- 3)ps {q - 2){q - 3){q - 4) {q - 3)p,, {q - 3)pio \ 

-{q-2){q-3) {q-2){q-3) 2(3 - g) 3-q 

q-2 (q-2){q-3) q''-7q + 13 3-q 

\-(q-2)(q-3) (g-2)(g-3) (q-3){q-5) (q - 3f J 



(5.14) 



where ps, etc. are defined in the appendix. This yields the same characteristic polynomial 
and eigenvalues as those obtained for the free Ly — A strip in Ref. [23] and can be taken as 
the transfer matrix for that strip. We have 
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det{Ttri,4,o) = (g - 2)%q - 3)' 



- lOg^ + 42g2 - 88g + 76 



(5.15) 
(5.16) 



(We recall that the actual form of the matrix is basis-dependent; for example, for the free 
Ly = 4 strip, a different matrix which, however, also has the same characteristic polynomial 
and eigenvalues originally calculated in Ref. [23], was given in Ref. [59].) 



tri,4,l — 





-^34 


-q2Pio 


92,3 


-92 


-93,5 


-92,3 


-Pl3 


-9i 


93 \ 




PlO 


-q2Pio 


-9293 


92,3 


93 


-92,3 


293 


-9i 


93 




-qs 


92,3 


-929I 


-92P1O 


93 


-92,3 


-Pl3 


-4 


-93,4 




1 


-92 


92,3 


-92P1O 


93 


-92,3 


-Pl3 


93 


-93,4 




qA 


92 








-93 
















-2 


292 


-92,3 


-92,3 


-93 


92,3 


-293 


-93 


-93 




-2 


92 


92 


92 


-293 


92,3 


Pis 


-93 


-93 




2^3 


-92,3 


-92,3 


292 


-293 


92,3 


Pu 


9i 


-293 













92 








1 





-93 / 



We calculate 



det{Ttri,4,i) 
Tr{Ttri,4,i) = 



(^_2)i2(^_3)8 
2(g-3)(2g2-12g + 21) 



triA,2 



/ 


PlO 


93,4 


-94 


92,3 


-92 





93 


\ 




-93 


93,4 


Pl4 


92,3 


92,4 


-92 


93 


93 




1 


-94 


Pu 


-92 


92,4 


-92 


93 


93 







-93 


-94 


92,3 


92,4 


92,3 





93 







1 


-94 


-92 


92,4 


92,3 





93 










1 





-92 


92,3 





93 




-2 


94 


94 


92 


92 





-93 





\ 








-2 





92 


92 





-93/ 



det{Ttri,4,2) = (9 - 2)^(9 - 3)^ 

Tr{Ttri,4,2) = Qq^ - 38g + 62 
As a special case of our general result in [39] , we have 

f q-3 q-A q-A q - 2\ 
—1 q — A q — A q 



-tri,4,3 — 



\ 
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- A q 
-1 q 



(5.17) 



(5.18) 
(5.19) 



(5.20) 



(5.21) 
(5.22) 



(5.23) 



The determinant and trace are the L„ 



2 
2 

-1 q-2j 
A special case of our general formulas (3.7) and 



(3.8), viz., det{Tt„,4,3) = (9 - 2)2(g - 3)^ and Tr{Ttri,4,3) = 13 - Aq. 
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J^. Ly = 5 



We present here new results for P{tri, 5 x m, BC, q) for BC=cyclic and Mobius. For 
d = we calculate 



tri,5,0 — 





q2S82 


g2,3Pl3 


92^^58 


Sl73 




92,3Pl3 


Sl48 


93^^34 


qsPn 


959i \ 




-q2Pw 




-292,3 


-9l 




-92,3 


-Pl3 


-93,4 


-93 


9i 






q2qi 


92P13 


-9l 




-92,3 


293 


93 


-93 


9i 




-q2 


-92,3 


92P13 






-92,3 


-Pl3 


93 


93,5 


ql 




-q2ql 


q2qi 


9294 


-94(29 - 


7) 


92 9I 


-H 


-9| 


93,4 


ql 




q2,3 


-92,3 


92P13 






929i 


qspis 


-9| 


93(29-7) 


ql 




-q2Pio 


-92,3 


-292,3 


-94(29 - 


7) 


929i 




93P10 


93,4 


ql 










-92 


-94 













93 





\ 


2^2 


92,3 


292,3 


94,5 




92,3 


Pl4 


-293 


293 


-93V 



(5.24) 

where the various p, r, and s polynomials are defined in the appendix. We find 

det{Ttri,5,o) = -(9 - 2)^^(9 - 3)^2 (5.25) 



Tr{Ttri,5,o) = (9 - 3) (9' - IO9' + 46g2 - lUq + 118) . (5.26) 

The characteristic polynomial of Ttri,5fl consists of factors of degree 2 and 7. 

The transfer matrices Ttri,5,d for li = 1,2 are both of dimension 21x21 and are too lengthy 
to present here; they are available from the authors. However, we do note the simple results 

detiTtri,,,i)^{q-2f%-3y' (5.27) 
Tr(Ttri,5,i) = -5?^ + 62?^ - 306?^ + 708? - 642 (5.28) 

det{Ttr^,,,2) = {q-2y'{q-3f^ (5.29) 

Tr{Ttri,5,2) = lOq^ - 98g^ + 332g - 387 . (5.30) 

Factorizations of CP{Ttri,5,j, z), j = 1, 2, are given in Table VIII. 
For d = 3 we calculate 
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tri,5,3 — 





-Pio 


-^3,4 


54 


-^3,4 


94 





-92,3 


92 








-93 





\ 






-?3,4 


-Pl4 


-?3,4 


-d 


94 


-92,3 


-92,4 


92 





-93 


-93 







-1 


Q4 


-Pl4 


Q4 


-d 


94 


92 


-92,4 


92 





-93 


-93 










13 


94 


-?3,4 


-94^ 


-Pu 


-92,3 


-92,4 


-92,4 


92 





-93 


-93 







-1 


?4 


?4 


-qI 




92 


-92,4 


-92,4 


92 





-93 


-93 










-1 







-Pu 





92 


-92,4 


92 





-93 


-93 













?3 


94 


94 


-92,3 


-92,4 


-92,4 


-92,3 








-93 













-1 


94 


94 


92 


-92,4 


-92,4 


-92,3 








-93 
















-1 


94 





92 


-92,4 


-92,3 








-93 



















-1 








92 


-92,3 








-93 




2 


-?4 


-?4 


-Q4 


-94 





-92 


-92 








93 
















2 





-94 


-94 





-92 


-92 








93 





I 

















2 








-92 


-92 








93 / 



(5.31) 

We calculate 

det{Ttr^,,,3) = (9 - 2)^^(9 - 3)^' (5.32) 



(5.33) 



The characteristic polynomial of Tj^i^s^s has factors of degree 5 and 8. 
For d — 4, our general theorem in [39] yields 



■ tri.bA 



/9 


-3 


g - 4 


9 


-4 


9 


-4 


9- 


2\ 




-1 


g - 4 


9 


-4 


9 


-4 


9- 


2 







-1 


9 


-4 


9 


-4 


9- 


2 












-1 


9 


-4 


9- 


2 

















-1 


9- 


2/ 



In accordance with our general results (3.7) and (3.8), we have 

det{Ttri,5A) = (9 - 2)'(9 - 3)=^ 

Tr{Ttri,,,4) = 5g - 17 . 
The characteristic polynomial of Ttri,5,4 has factors of degree 2 and 3. 



(5.34) 



(5.35) 
(5.36) 
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VI. STRIPS OF THE HONEYCOMB LATTICE 



A. Ly = 2 

The chromatic polynomials P{hc, 2 x m, BC, q) for BC=cychc and Mobius were given in 
[27]. For d = 0, Thc,2fl is a scalar, Thc,2fl = Kc,2fi ^ Dq ^ q"^ - 5q^ + lOq^ - lOq + 5. For 
d — 1, from our general formulas above, we obtain 

2-q 1 
1 2-q \ (6.1) 
-1 -1 



1-q -1 

^hc,2,i,2 =1 1 - g -1 I (6.2) 
1 



whence 

_/D4 -DA _ f q'-3q + 3 2-q \ 
^"'''^'-[-D, D, )-[ 2-q q^-3q + 3) ^^'^^ 

yielding the relevant special case of the eigenvalues given in eq. (6), (7) of [27], Xhc,2,i,i — 
Di-D2.^ ^4,3 ^q^ -Aq + b and \hc,2,i,2 = -D4 + -D3 = (g - 1)^. We note that 



det{Thc,2,i) = {q- 1)1^ - 4g + 5) (6.4) 



and Tr{Thc,2,i) = 2L>4 



B. Ly = 3 



Further, 



so that 



qi,2S7 ?2S3,5 i5,5 

Thc,3,o- I ?i ^4,3 -qs I (6.5) 

--^"6,2 -?2P4 -i^5,4 



c?ei(r^c,3,o) = (g-i)'(g-2)2 (6.6) 

Tr(nc,3,o) = - + 28g^ - 56q^ + 71q^ - 58q + 26 . (6.7) 



34 



For d = 1, 



T; 



hc,3,l 



( ?1,2P4 -?l-^4,3 1 ''IS -^4,3 ?2P6 \ 

gi,2^4 -^5 -^2,3 ^3-04 -2g| 

91,2 -giA ^6,5 -92,3 93^ 92^4,3 

-gi,2 qi -gg 1 -g2 

-qi,2 qiDi Db -2^2 2^4 52,3 

V gi qi -1 -?3 -qs -^4,3 / 



from which it follows that 



(6.8) 



det{Thc,3,i) = {q-lT{q-2f 



(6.9) 



Tr{Thc,3,i) = - W + 46g2 - 60g + 37 . 

For d — 2, applying our general formulas above, we calculate 

/2-q 1 \ 

1 2-q 0-1 

^^.0,3,2,1 = 1 - g -1 

-1-1 

\ 1 y 



(6.10) 



(6.11) 



S. 



hc,3,2,2 



(l~q -1 0\ 
2-q 1 -10 
1 2-q 
1 
\ -1-10 0/ 



(6.12) 



The factorizations for A = he arc given in Table IX. We note that for even Ly, the factor- 
izations in this table correspond to the numbers np{hc, Ly, d, +) and np{hc, Ly, d, —) given 
in Table 11, when these numbers can be defined, but there are no nontrivial factorizations 



for odd Ly. We calculate 



det{Thc,3,2) -{q- 1)' 



(6.13) 



Tr{nc,3,2)=3q'-I0q + 12 



(6.14) 
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C. Ly=4 

For the 6x6 dimensional matrix Thc,4,Q we calculate 

det{Thc,Afl) = 2(g - lf{q - 2f{q^ - 5g + l){2q^ - 20q^ + 83?^ - 185?^ + 23 V - 175g + 60) 

(6.15) 

Tr{Thc,A,o) =q^ -iW + 55q^ - 165q^ + 333g^ - 480g^ + 503^^ - 362g + 142 . (6.16) 

The characteristic polynomial has factors of degree 1 and 5, as indicated in Table IX. The 
eigenvalues are 

A;,c,4,o,i = (?-l)'(g'-5g + 7) (6.17) 

and the roots of the factor of degree 5. 

For the 13 x 13 dimensional matrix 7)10,4,1 we calculate 

det(Thc,4,i) = 4(g - l^^q - 2f{q^ - Ibq^ + lOOg^ - 389g^ + 974g^ - 1623g^ 

+1774g2 - 1171g + 360) (g^^ - 20g^^ + 185g^° - 1048g^ + 4061g^ 

-11385g^ + 23784g6 - 37468g^ + 44352g^ - 38630?^ + 23631?^ - 9200g + 1750) 

(6.18) 

Tr{Thc,4,i) = 4g^ - 36g^ + 147g^ - 351?^ + 529?^ - 489g + 229 . (6.19) 

The characteristic polynomial factorizcs into factors of degree 5 and 8. 
For the 11 x 11 dimensional matrix Thc,A,2 we find 

det(Thc,4,2) = (g - iy^(q^ - lOq^ + 44?^ - 108g^ + 160?^ - 141g + 60) 

X (2g^ - 29g^ + 186g^ - 690g^ + 1634?^ - 2564g=^ + 2647?^ - 1670? + 500) (6.20) 

Tr{Thc,4,2) = Gq^ - 39g^ + llOg^ - 157? + 104 (6.21) 
For Thc,A,3, applying our general formula (3.12), we calculate 

36 



T, 



/ic,4,3 



whence 



The eigenvalues are 



f D4 2-q 0\ 

3-g F4,3 2-g 1 g-3 
1 2 - g F4,3 3 - g g - 3 

2 -g D4 
\-lq-2q-2-l 2 / 



det{nc,4,3)^2(q-l)\q^-Aq + 5) 



rr(7),e,4,3) = 2(2g2-7? + 9) 



A/ic,4,3,l — (9 - 1)' 



A 



hc,4,3,2 



Fa 3 = g - 4g + 5 



and the three roots of 



- 2(g2 - 4g + 6)^^ + (g^ - + 26g^ - 36g + 21)z - 2{q - 1)^ = 



D. = 5 



We have calculated the T^^^^^^ and find 

det{nc,^,,)^{q-lf\q-2f^ 

Tr(Thc,5,o) = - 1 V + 91g^ - 364g^ + 1007?^ - 2058?^ + 3232g^ 
-3950g^ + 3662^2 - 2350g + 808 



det{Thc,5,i) = (g - lY'iq - 2) 



\24 



Tr{Thc,5,i) = 5q^ - 60g^ + 336g^ - 1152g^ + 2676g^ - 4374g2 



+499 V - 3720? + 1429 
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\16 



det{THc,5,2)^{q-ir{q-2) 



Tr{{Thc,5,2) = lOg^ - 96q^ + 422g^ - 108 V + 1742^2 - 1686g + 805 



det{n,,,,s)^{q-iy"{q-2) 



Tr{Thc,5,3) = lOg^ - 6Sq^ + 202g=^ - 302g + 207 



(6.32) 
(6.33) 
(6.34) 
(6.35) 



det{Tnc,5A) = iQ- 1)' 

Tr(T;,e,5,4) =5g2-18g + 24 
The matrix Thc,5,4. is sufficiently small that we can list it here: 



/L'4 2-g0 0\ 

3-q Fi^3 2-q 1 q-3 

1 2-q F4,3 3-q q-3 

3 - ? F4,3 1 - ? q-3 

1 2-q gi,2 g-3 

-1 g-2 g-2 -1 2 

\0 -lq-2q-lQ 2 ) 



hc,5,4 — 



(6.36) 
(6.37) 



(6.38) 



Factorizations of characteristic polynomials are given in Table IX. 



VII. SOME ALGEBRAIC PROPERTIES OF EVALUATIONS OF TRANSFER 

MATRICES 

We have studied algebraic properties of evaluations of transfer matrices T\^Ly,d at special 
values of q. We find many interesting results, and mention only a few here. 

First, from our exact results for the Ta^Lj,,^ for c? = Lj^ — 1, we can calculate properties of 
powers of these matrices, in particular, when evaluated at special values of q. For example. 



m r)m— 1 ' 1 1 

sg,2,lj - ^ \ I I 



for q — 1, n > 1 



(7.1) 



and, for odd m = 2n + 1 and even m — 2n, 



1 

1 



[i;,,2,i]'"+^ = ( ; ; ) for g = 2, n > 



(7.2) 
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[T.,,2,i]'" = ( J ? ) ? = 2, n > 1 (7.3) 

and 

[Ts,,2,ir = (-l)'"2™-i Y ) ^ = 3' ^ > 1 ■ (7-4) 

Two other identities are given below (where n > 1 and there is no sum on j): 

m,,,,,n, = 2^-\l + 2--') for g = 3 (7.5) 

[(T.,,6,5)n-, = ^ [2^--^ + 3" + l] for g = 3. (7.6) 

We next consider transfer matrices for strips of the triangular lattice. In general, from eq. 
(3.7), it follows that the evaluation ofTtri,Ly,Ly-i at g = 2 is noninvertible and the evaluation 
of this matrix at g = 3 is noninvertible if L^, > 3. First, we observe that 

[Ttri,2,ir = ( 1 ) for g = 2 . (7.7) 
The matrix Ttri^2,i at g = 3 and its powers form the cyclic group of order 3, Z3, with 

[Ttri,2,ir+^ = [Ttri,2,ir ^nd 

T -f^-^\ IT IT 13-/^10 

J-tri,2,l ~ \ I _l j ' [■'■tri,2,l\ - [ _i q j > L^tH,2,lJ " I Q 1 

for q = 3 (7.8) 

Given that the transfer matrices are real, if such a matrix satisfies (TajL^.^)^ = Ik, where 1^ 
denotes the k x k identity matrix, this implies that det(Tjy Ly,d) — il) ^-nd if P is odd, as is 
the case here, with p — 3, then det(T\^L^ (i) = 1- Referring to our general result (3.7), we see 
that this condition is, indeed, met with A; = 2 if L^^ = 2, = 1, and g = 3. 

The matrix Ttri,3fl evaluated at g = 2 satisfies the relation (denoting this by T for 
short) = T, so, together with the 2x2 identity matrix I2, {Ttri.3fl)q=2, [{Ttri,3,o)q=2]'^ 
form a multiplicative semigroup of order 3. (Here, we recall that a semigroup satisfies the 
same axioms as a group but without the axiom that inverses of elements must exist; this is 
also called a semigroup with unit or monoid [67]. The same matrix evaluated at g = 3 is 
idempotent, i.e., satisfies = T, so that I2 and (Ttri,3,o)g=3 form a semigroup of order 2. One 
surmises that connections can be made between these cycle lengths and the dependence of 
the chromatic number of the cyclic strip of the triangular lattice on its length (see appendix) , 
but we have not pursued this. 
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The matrix Ttri,3,i evaluated at g = 3 satisfies the relation (again denoting this by T 
for short) T"^ = T so that h, {Ttri,3,i)q=3, [(^^r^,3,l)<7=3]^ and [(Tt^i,3,i)g=3]^ form a semigroup 
of order 4. The matrix Ttri,3,2 evaluated at g = 3 satisfies the relation = T, and hence 
the elements I3, {Ttri,3,2)q=3, and [{Ttri,3.2)q=3]'^ form a semigroup of order 3. (These are not 
groups since {Ttri,3,i)q=3 and {Ttri,3,2)q=3 are not invertible.) We have found a number of 
similar algebraic properties of powers of the T\^Ly,d matrices evaluated for special values of 

VIII. DISCUSSION OF DETERMINANTS 
A. Square-Lattice Strips 

Here we discuss further the determinants det{T\^Lyj) for A = sq.tri. hc, starting with 
the square lattice . We first consider the zeros of the determinants, restricting to the range 
< d < Ly — 1 since Tsq^Ly,Ly = (—1)^^ has no zeros. From our theorems in Ref. [39], it 
follows that det{Tsq^Ly,Ly~i) has (i) the factor (g — 1) always; (ii) the factor (g — 3) if L^, = 
mod 2; (iii) the factors (g — 2) and (g — 4) if = mod 3; further, all of the zeros of 
det{Tsq^Ly,Ly-i) lie in the interval 1 < g < 5 and, as Ly 00, they become dense in this 
interval. One may also investigate the properties of the zeros of Tsq^i^^a ioi other values of 
d than Ly — 1. Below we list some zeros; for multiple zeros, we indicate the multiplicity in 
parentheses, (m = 2), etc. 

det{Tsq,2,o) = at g = -(3 ± iv^ ) = 1.5 ± 0.866i (8.1) 
det{Tsq,3,o) = at g = 1, 1.659 ± 1.1615i, 2.682 (8.2) 
det{Tsq,3,i) = at g = 1, 1.229 ± 0.712i, 2 (m = 2), 3.271 ± 0.449i (8.3) 

det{Tsq,3,2) = at g = 1, 2, 4 (8.4) 

det{Tsq,ifi) ^ Q at g = 1, 1.163 ± 0.539i, 1.355, 1.730 ± 1.331i, 2.437, 3, 3.210 ± 0.6918i 

(8.5) 

det{Tsq,4,i) = at g = 1 (m = 3), 1.257 ± 0.578i, 1.302, 1.5 ± 0.866i (m = 2), 
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1.807 ±0.936i, 2.736, 3(m = 5), 3.4169 ± 0.6339i 



(8.6) 



det{Tsg,4,2) = at ? = 1 (m = 2), 1.180 ± 0.600i, 1.327, 1.580 ± 0.606i, 



2.198 ±0.573i, 3(m = 2), 3.839, 3.959 ±0.294i 



(8.7) 



det{T,g^4,3) = at q ^ 1, 1.586, 3, 4.414 



(8.8) 



and so forth for higher Ly values. We plot these zeros in Figs. 4-5. (The figures are formatted 
such that, to save space, the vertical axis is shown intersecting the horizontal axis slightly 
to the left of g = 1.) 

We observe that all of the zeros of det{Tsg Ly fi) that we have calculated satisfy the con- 
dition Iq' — 3| < 2. As is evident in eq. (8.15), some zeros of det{Tfif.Lyj) have real parts 
that lie outside the interval 1 < Re{q) < 5. It would be of interest to prove a theorem 
bounding the complex zeros of the determinants analogous to results that have been ob- 
tained for chromatic zeros, such as the theorem in Ref. [68]. We also note that for c? 7^ 
(as well as d ^ L), for the cases that we have calculated, det{Tsq,Ly,d) has the factor (g — 1). 
Although detiTsqfifi) = Tsq^^ does not have this factor, we observe that for the cases we 
have calculated with Ly > 3, det{Tsq^Ly,o) contains the factor (g — 1). 



2 




lm(q) 1 







-1 



-2 



FIG. 4. Zeros of det{Tsq^L) for L = 4 in the complex q plane. 
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FIG. 5. Zeros of det{Tsq^L) for L = 5 in the complex q plane. 

B. Triangular-Lattice Strips 

An interesting finding is that the determinants of Ttri,Ly,d are the simplest among the 
regular homopolygonal lattices, square, triangular, and honeycomb. Restricting to the range 
< d < Ly — 1 to exclude Ttri,Ly,Ly = (^1)^^, which has no zeros, we observe that all of the 
results that we have obtained are consistent with the following formula, which we conjecture 
to hold for general Ly (over the range Ly >2 where the triangular lattice strip is defined): 

det{Ttri,Ly,d)-ri{q-mQ-^Y (8-9) 

where r] = ±1 and 

a = 2np{tri, ,Ly-l,d) + {Ly - 2)np{tri, Ly - 2, d) (8.10) 

and 

b= {Ly-2)np{tri,Ly-l,d) . (8.11) 

Here the sign rj depends on Ly and d, and we take np{tri, l,d) = np{sq, l,d), in accord 
with the general equality np{tri, Ly,d) — np{sq, Ly, d) [9]. (It is not necessary to give a 
formal definition of np{tri, 0, d) because the only time it appears in the equation for a, for 



42 



the lowest nontrivial case Ly = 2, it is multiplied by zero.) Some structural properties of this 
conjecture are indirectly motivated by the exact formula that we have found in Ref. [66] for 
det{Tz^tri,Ly,d), where Tz,tri,Ly,d is the transfer matrix, in the degree-rf subspace, for the full 
Potts model partition function on the given strip. However, aside from the factors of {v + 1) 
in Tz^tri,Ly,d that cause it to vanish in the special case v = —1 corresponding to the chromatic 
polynomial, the rest of the expression that we have established for det{Tz,tri,Ly,d) does not 
involve factors that reduce to powers of (g — 2) or (g — 3) when v = —1\ instead, it involves 
powers of (g + v) and f , which reduce to (g — 1) and ±1. The main structural feature of our 
known result for det{Tz,tri,Ly,d) that influenced our present conjecture is that the powers of 
these factors involve quantities such as nz{Ly,d) and nz{Ly — l,d), where nz{Ly,d) is the 
dimension of Tz,K,Ly,d (which is the same for A = sg, tri, he) [9,10]. This led us to construct 
and test conjectured forms for the chromatic polynomial det{Ttri,Ly,d) that involve powers of 
the analogous dimensions np{tri, Ly,d). 

The conjecture (8.9) imphes that the zeros of these determinants det{Ttri,Ly,d) occur at 
the two points g = 3 and g = 2 (for < d < L^^ — 1), in accordance with all of our explicit 
calculations. It is interesting to observe that these are, respectively, the chromatic number 
xitri) and x(^^0 — 1 of the two-dimensional lattice. For the special case d — Ly — 1, eq. 
(8.9) (with the plus sign) reduces to our previously proved result (3.7). For Ly > 2 and 
d — Ly — 2, eq. (8.9) can be written explicitly as 



Starting from the other end of the range of d, for d = 0,1 (and Ly > 2), eq. (8.9) can be 
written explicitly in terms of the Motzkin numbers Mn defined in eq. (2.12) as 



det{T,„,Ly,Ly-2) = V{q - 2)^^^-^(g - 3)(^«-^)(^^-2) 



(8.12) 



det{Ttri,Ly,o 



(g - 2)^^^!'-2 + (^2'~2)^^!/-3(g - 3)(^!/-2)A^Lj,-2 



(8.13) 



det{Ttri,Ly,i 



) =n{q- 2)2^^y-i+(^^'~2)^^!/-2(g - 3)(^!/-2)Mij,-i 



(8.14) 



Similar explicit formulas can be obtained from (8.9) for other values of d. 



C. Honeycomb-Lattice Strips 



In the case of the honeycomb lattice, we have 



det{Thc,2,Q) ^ Q at 0.691 ± 0.951i, 1.809 ± 0.588i 



(8.15) 



det(Thc,2,i) =0 at 1 (m = 2), 2 ± i 



(8.16) 
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As is evident in our results above, for < < 2, the zeros of det{Thc,3,d) occur at the two 
points q = 1 and q = 2. Restricting to < < — 1 to exclude the constant Thc,Ly,Ly = li 
and to odd widths Ly, we observe that all of our results are consistent with the following 
formula, which we conjecture to hold in general: 



where 



and 



det{Thc,Ly,d)-iQ-'^yiQ-'^f for odd Ly (8.17) 



{Ly + l)np{hc,Ly-l,d) (8.18) 



t= {Ly-l)np{hc,Ly-2,d) (8.19) 



where we take np{hc,l,d) = np{sq,l,d). The motivations for this conjecture are based 
in part on the exact expression that we have found elsewhere [66] for det{Tzjic,Ly,d), where 
Tz,hc,Ly,d is the transfer matrix, in the degrec-d subspace, for the partition function of the full 
Potts model on the given strip. However, we note that our exact expression for det(Tz^hc,Ly,d) 
docs not have a factor that reduces to a power of (g — 2) for the special case = — 1 that 
defines the chromatic polynomial. Thus, as in the case of the other lattices, the detailed 
structure of the determinants of Tz,K,Ly,d for the full partition function and Tp^/<^^Ly,d, denoted 
here simply as T^^Ly^d, for the special case v = —1 are different. This is not surprising, 
since the matrix Tp_A,Lj,,d does not arise, in general, from a block decomposition of Tz^A,Ly,d 
but instead by removing zero columns and corresponding rows of the latter matrix. It 
is interesting that the zeros of det{Thc,Ly,d) occur at q = x{hc) and q = x{hc) — 1, where 
x{hc) = 2 is the chromatic number for the honeycomb lattice. For the special case d — Ly — 1, 
eq. (8.17) agrees with the result (3.30) that we have proved from our general calculation of 
Thc,Ly,Ly-i- For the special case d — Ly — 2, eq. (8.17) can be written exphcitly as 

det{Tnc,Ly,Ly-2) = {q-lf''^^'^^''''-'^'\q-2)'^y-' for odd Ly . (8.20) 

Similar exphcit formulas can be obtained from (8.17) for other values of d. Both of our 
conjectured general formulas for det{Tp^^Ly,d) have the form 



pi 



P2 



det{TA,Ly,d)-r][Q-{x{^)-^ - x(A)J for A = iri,/ic (8.21) 

where x(A) is the chromatic number of the respective two-dimensional lattices A — tri, he; 
77 = ±1 for A = tri and 77 = 1 for A = /ic; Lj^ > 3 is odd for the honeycomb-lattice strips; 
and the powers pi and p2 have been given in (8.10), (8.11), (8.18), and (8.19). 
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IX. ACCUMULATION LOCUS OF CHROMATIC ZEROS 



Although the locus B for a strip with a given width Ly and free longitudinal boundary 
conditions (and free or periodic transverse boundary conditions) is different from that for the 
corresponding strip with periodic or twisted periodic longitudinal boundary conditions, one 
notices that in some ways these loci become more similar as Ly increases. One expects, for 
example, that in the limit Ly — > oo, with any of these boundary conditions, one will obtain 
the same value of qc for the B for the Potts model on the resultant infinite 2D lattice. Prom 
our new exact calculations, we find the following results. Pirst, for the 5 x oo cyclic/ Mobius 
strip of the triangular lattice, B crosses the real q axis at q — 0,2,3, and the maximal point 

qc = 3.33245.. , A = tri, Ly x L^ ^ 5 x oo, BC = cyclic/Mobius . (9.1) 

Prom our new calculations for the Ly — A and Ly — 5 cychc/ Mobius strips of the honeycomb 
lattice and our analysis of the locus B in the Lj; — > oo limit, we find that in both cases B 
crosses the real q axis at the points q — and q — 2, and at the maximal crossing points 

qc = 2.15476.. , A ^ he, Ly x L^ ^ A x oo, BC = cyclic/Mobius (9.2) 

qc = 2.26407... , A^hc, LyX L^^bxoo, BC = cyclic/Mobius . (9.3) 

In Table X we present a summary of these new results, together with the results that we have 
obtained for smaller values of Ly in earlier work. Our present calculations are in agreement 
with our previous conjecture, based on all of the exact results that we had obtained, that 
for the infinite-length limit of a cyclic or Mobius strip graph of a given lattice A, with a 
given width Ly, qc is a nondecreasing function of Ly. For each type of lattice, the values of 
qc increase toward the values for the infinite 2D lattices, qdsq) = 3 [40], qdtri) = 4 [42], 
and, formally, qdjic) = (3 + \fE)/2 ~ 2.618 [48,43]. Although the values of qc for each of 
these types of lattice strips are less than the value of qc for the respective two-dimensional 
lattices, we have shown previously that this need not necessarily be the case. Por example, 
for (the — > oo limit of ) the Ly = 3 strip of the square lattice with toroidal boundary 
conditions [32], qc — 3, equal to the value for the square lattice, and, indeed, for the cyclic 
self-dual strip of the square lattice, ?c = 3 for each of the widths for which we obtained exact 
calculations [69,11,70]. Nevertheless, the observed monotonically nondecreasing behavior of 
the qc values provides a nice interpolation between the value qc — 2 for the one-dimensional 
Potts antiferromagnet and the values on the respective two-dimensional lattices. As noted 
before, in general, if one uses strips with free longitudinal boundary conditions, the locus B 
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does not necessarily cross the real q axis [23], and in cases where it does not, there is, strictly 
speaking, no qc{{G}). In these cases, the way that one gains information is different; one 
examines the endpoints of the rightmost arcs on B that are closest to the real axis. These 
are found to exhibit the general tendency to move gradually to the right with endpoints 
that move closer to the real axis as the strip width is increased [23,24,34,57-59,35,39,36]. 
The results are consistent with the inference that as Ly oo, the arcs will merge to form 
closed curves on B and the right-most part of B will cross the real axis at the value of Qc for 
the corresponding two-dimensional lattice. (One can attempt to carry out a similar type of 
study for three-dimensional lattices, and exploratory work has been done in [38].) 

In Figs. 6-9 we show the singular loci B for infinite-length strips of the triangular lattice 
with width Ly — 5 and of the honeycomb lattice with widths Ly — 3, 4, 5, together with zeros 
of the chromatic polynomials (chromatic zeros) for long finite-length strips of each type. We 
first discuss the triangular- lattice strip (Fig. 6). This should be compared with our earher 
results for the Ly — 2 strip in Fig. 3 in Ref. [30] and for the Ly — 3 and Ly — A strips in 
Figs. 2 and 3 of Ref. [36] . As expected, the comparison is particularly close with the Ly — A 
case. For L^^ = 3, 4 we found that there are two inner curves that run between points on the 
outer envelope and cross the real q axis at g = 2 and q — S, separating the interior of the 
envelope curve into regions that include the real segments (i) < g < 2, (ii) 2 < g < 3, and 
(iii) 3 < q < qc- We also found that the densities of zeros on these curves passing through 
q — 2 and q = 3 were slightly less than the densities of the zeros on the outer envelope 
curves. Another feature that we found was the existence of two "bubble" regions protruding 
to the right from the main curve that passes through qc- Finally, we observed that as Ly 
increases from Ly = 3 to Ly = 4, there are more zeros, and more support for a part of B in 
the half-plane with Re{q) < 0. Our current results for Ly = 5 exhibit these same features, 
and show an increase, relative to our Ly = 4 results, in the portion of B extending into the 
half-plane with Re{q) < 0. This is indicated in Tabic X. In this table we also characterize 
the form of the curve B as it passes through the point q = 0; given that it passes vertically 
through this point, and is invariant under complex conjugation, there are two possibilities, 
viz., that it extends into the half-plane with (i) Re{q) > or (ii) Re{q) < 0. We denote these 
two possibilities as (i) convex and (ii) concave to the left. Note that, a priori, the conditions 
that B has support in the half- plane with Re{q) < and that it is concave to the left are 
not equivalent since it might be convex to the left, i.e., curve into the Re{q) > half-plane 
in the neighborhood of the origin, but then curve back to the left at larger values of \Im{q)\ 
and cross over into the Re{q) < half-plane. However, we find that this does not happen 
for any of the strips for which we have done calculations, as is evident in Table X. Thus, for 
these strips, the properties that B has support for Re{q) < and that it is concave to the 
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left at g = are observed to occur together. We also note that the zeros on the right-hand 
side of the locus B cross the real axis at a point consistent with the asymptotic result for 
given in eq. (9.1). We observe small complex-conjugate bubble regions around q = 2.6 ± 2i, 
where the curve passing through q = 2 intersects the outer curve on B. We should also 
remark on the parts of the locus B near the points q = 1 ±2.1i: without magnification, these 
appear to be arc endpoints, but in fact the end in small, very narrow, bubble regions. This is 
consistent with our finding for the ^ oo limit of all strips of regular lattices with periodic 
longitudinal boundary conditions, that the respective loci B do not have arc endpoints, in 
contrast with the situation for strips with free longitudinal boundary conditions, for which 
the loci B generically do exhibit such arc endpoints. In our earlier work we have found 
that loci B can exhibit tiny sliver regions, and we mention that for this case of the Ly — 5 
cychc strip of the triangular lattice and the for the loci shown below for the strips of the 
honeycomb lattice, such tiny sliver regions, if small enough, would elude our analysis because 
of the finite size of the grid that we use for testing for equimodular dominant A's. 

In addition to the continuous accumulation set of chromatic zeros that forms B, one may 
also investigate discrete chromatic zeros. In general, for any graph with at least one edge, 
there is always a chromatic zero at q — 1, since there is no proper coloring of such a graph 
with just one color. For strips of the triangular lattice with width the Ly = 5 and reasonably 
great length, we also find a chromatic zero near to the point g = (3 V5)/2 ~ 2.618, which 
is a zero of c^'^\ where a degree c? = 2 A is dominant, and a chromatic zero near to the point 
q ~ 3.2470, a zero of c'^^\ where a degree = 3 A is dominant. 

We next show the singular loci B for the infinite-length strips of the honeycomb lattice 
with Ly = 3,4,5 in Figs. 7-9, together with chromatic zeros for long finite-length strips of 
each respective width. As background, we recall that in the ^ oo limit for the smaller 
widths Ly = 2 (Ref. [27]) and Ly = 3 (Rcf. [10]), wc found that six curves on B, forming three 
branches, intersect at the point qc (which is equal to 2 for those widths). Hence, it was found 
that the region diagram included at least six regions (which comprised the totality of regions 
for Ly — 2 and were augmented by two additional very small regions centered approximately 
at q — 0.5±0.45i for Ly — 3): (i) the outermost region Ri, extending infinitely far away from 
the origin and including the intervals q > 2 and q < on the real axis; (ii) the innermost 
region i?2, which includes the interval < q < 2; (iii) a complex-conjugate (c.c.) pair of 
regions R3, forming upper and lower outer crescent-shaped areas adjacent to qc, and (iv) 
the c.c. pair R^, Rl forming upper and lower inner crescent-shaped areas adjacent to qc- It 
was also found that the density of zeros on the inner two of the six curves passing through 
qc was somewhat lower than the density of the zeros on other parts of B. The locus B in the 
vicinity of the origin q — is concave to the left, which implies that it has support for values 
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of q in the half-plane with Re{q) < 0. This concavity increases as the width increases. In 
the loci B presented in Fig. 8 and 9, we see that the point g = 2 is an intersection point of 
four rather than six curves, and this feature is correlated with the fact that for these widths, 
Qc exceeds 2. Thus, in the plot for = 4 we see a small self-conjugate bubble region that 
includes the interval 2 < g ^ qhc4.c, where qhc4.c — 2.155 is the value of Qc for this strip, given 
in eq. (9.2). For Ly = 5 we observe three small self-conjugate regions in B (sec Fig. 9) that 
include the three segments 2 < q < qhc5a, qhc5a < q < qiu^b, and qhc^b < q < qhc5c, where 
qhc5a — 2.1997, qhc5b — 2.2468, and qhcbc — 2.26407 is the value of qc for this strip, given in 
eq. (9.3). Thus, the locus B for the cyclic Ly — 5 strip of the honeycomb lattice crosses the 
real q axis q = 0, 2, qhc5a, Qhcbb, and qhcbc- These values are listed, together with those for 
other strips, in Table X. For the Ly — 3 strip, one sees in Fig. 7 a complex-conjugate pair 
of narrow regions around g ~ 0.5 ± lAi; similarly, for the Ly — A strip, one sees in Fig. 8 a 
complex-conjugate pair of very narrow sliver regions around g ~ 0.1 ± 1.3i. For the Ly — 5 
strip of the honeycomb lattice, we find arcs ending in tiny bubble regions at approximately 
g = —0.1 ± i. A general feature that we observe is increasing complexity of the loci B with 
increasing strip width. As before, we note that regions that are extremely small could be 
missed by our grid used in testing for degeneracies between dominant eigenvalues that yield 
components of the locus B. 
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FIG. 6. Singular locus B in the complex q plane for the ^ oo limit of the cyclic strip of the triangular 
lattice with width Ly = 5. For comparison, the plot also shows zeros of P{tri,Ly x Lx,cyc.,q) for Ly = 5 
and a typical large value of the length, = 20. 
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FIG. 7. Singular locus B in the complex q plane for the Lx ^ co limit of the cyclic strip of the honeycomb 
lattice with width Ly = 3. For comparison, the plot also shows zeros of P{tri, Ly x L^, eye, q) with Ly = 3 
and a typical large value of the length, = 40. 
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FIG. 8. Singular locus B in the complex q plane for the ^ oo limit of the cyclic strip of the honeycomb 
lattice with width Ly = 4. For comparison, the plot also shows zeros of P{tri, Ly x L^, eye, q) with Ly = A 
and a typical large value of the length, = 24. 
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FIG. 9. Singular locus B in the complex q plane for the ^ oo limit of the cyclic strip of the honeycomb 
lattice with width Ly = 5. For comparison, the plot also shows zeros of P{tri, Ly x L^, eye, q) with Ly = 5 
and a typical large value of the length, = 20. 
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X. APPENDIX 1: CHROMATIC NUMBERS FOR LATTICE STRIPS 



The following chromatic numbers apply for lengths that are greater than the lowest few 
where strips sometimes degenerate. For cyclic strips of the square lattice, the chromatic 
number is given by 

IT T \ f2 for even La; /-.^^ -,\ 

xN,L,xL.,c^c.) = |3 ^^^^^^^^ (10.1) 

independent of Ly. For the Mobius strips of the square lattice (e.g., [35]) 

^ ^ \3 for (La;, Lj,) = (e,e), (o,o) ^ > 

where e and o denote even and odd. 

The cyclic and Mobius strips of the triangular lattice have chromatic numbers (e.g. [36]) 

X(trz, Ly X La., eye.) = { ^ \ "^f^^^ ^ 3 (10.3) 

X{trh Ly X La;, Mh) = 4 (10.4) 

Again for L^ above the first few values [L^ = Ij 2 where strips can degenerate), the cyclic 
and Mobius strips of the honeycomb lattice have chromatic numbers (e.g. [10]) 

x{hc,Ly X L^,cyc.) ^2 (10.5) 
X{hc, Ly X La;, Mb.) = 3 (10.6) 
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XI. APPENDIX 2: TABLES 



TABLE I. Tabic of mimbcrs np{A,Ly,d) and their sums, Np^A,Ly,\ for cyclic strips of the square and 



triangular lattices. Blank entries are zero. 



Ly i d ^ 





1 


2 


3 


4 


5 


6 




1 


1 


1 












2 


2 


1 


2 


1 










4 


O 


2 


-1 


3 


1 








10 


4 


4 


9 


8 


4 


1 






26 


5 


9 


21 


21 


13 


5 


1 




70 



TABLE IL Table of numbers np{hc, Ly, d) and their sums, Np^hc,Ly,\ for cychc strips of the honeycomb 
lattice. Blank entries are zero. 



Ly i d ^ 





1 


2 


3 


4 


5 


6 


Np^hC,Ly,X 


2 


1 


2 


1 










4 


3 


3 


6 


4 


1 








14 


4 


6 


13 


11 


5 


1 






36 


5 


19 


43 


40 


22 


7 


1 




132 



TABLE III. Table of Anp{sq,Ly,d) for strips of the square lattice. Blank entries are zero. The last 
entry for each value of Ly is the total number of partitions with self-reflection symmetry. 



Ly i d^ 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




1 


1 


1 




















2 


2 


1 





1 


















2 


3 


2 


2 


1 


1 
















6 


4 


2 


1 


2 





1 














6 


5 


5 


5 


3 


3 


1 


1 












18 


6 


5 


3 


5 


1 


3 





1 










18 


7 


13 


13 


9 


9 


4 


4 


1 


1 








54 


8 


13 


9 


13 


4 


9 


1 


4 





1 






54 


9 


35 


35 


26 


26 


14 


14 


5 


5 


1 


1 




162 


10 


35 


26 


35 


14 


26 


5 


14 


1 


5 





1 


162 
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TABLE IV. Tabic of numbers np{sq, Ly,d,zt) for strips of the square lattice. For each Ly value, the 
entries in the first and second lines are np{sq, Ly, d, +) and np{sq, Ly, d, —), respectively. Blank entries are 
zero. The last entry for each value of Ly is the total Np^l^^\. 



Ly 


n 1 

u,+ 


1 1 

1, + 


2,+ 


3, + 


A 1 


1 


6,+ 


1 1 

7,+ 


Q 1 

8,+ 


O 1 

9, + 


1 n 1 






0,- 


1, - 


2,- 


3, - 


A 

4, - 


5, - 


6, - 


7, - 


8, - 


9, - 


10,- 


AT" 


2 


1 


1 


1 
























1 




















4 


3 


2 


3 


2 


1 






















1 


1 


















10 


4 


3 


5 


5 


2 


1 


















1 


4 


3 


2 
















26 


5 


7 


13 


12 


8 


3 


1 
















2 


8 


9 


5 


2 














70 


6 


13 


27 


30 


20 


11 


3 


1 














8 


24 


25 


19 


8 


3 












192 


7 


32 


70 


77 


61 


34 


15 


4 


1 












19 


57 


68 


52 


30 


11 


3 










534 


8 


70 


166 


199 


163 


106 


49 


19 


4 


1 










57 


157 


186 


159 


97 


48 


15 


4 








1500 


9 


179 


435 


528 


468 


318 


174 


72 


24 


5 


1 








144 


400 


502 


442 


304 


160 


67 


19 


4 






4246 


10 


435 


1107 


1405 


1288 


946 


550 


265 


96 


29 


5 


1 






400 


1081 


1307 


1274 


920 


545 


251 


95 


24 


5 




12092 



TABLE V. Table of Anp{hc, Ly, d) for strips of the honeycomb lattice with even Ly. Blank entries are 
zero. The last entry for each value of Ly is the total number of partitions with self-reflection symmetry. 



Ly i 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


^Np,hC,Ly 


2 


1 





1 


















2 


4 


4 


3 


3 


1 


1 














12 


6 


7 


-1 


7 


1 


1 





1 










21 


8 


36 


30 


30 


17 


17 


6 


6 


1 


1 






144 


10 


66 


47 


66 


23 


47 


7 


23 


1 


7 





1 


288 
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TABLE VI. Tabic of numbers np{hc, Ly,d,zt) for strips of the honeycomb lattice with even Ly. For 
each Ly value, the entries in the first and second lines are np{hc, Ly, d, +) and np{hc, Ly, d, —), respectively. 
Blank entries are zero. The last entry for each value of Ly is the total Np^hc,Ly,\- 



T ( rl 

Ly [d, +) 


n 1 
U, + 


1 1 


O 1 


Q 1 


A _L 

4, + 


K 1 

0,+ 


a 1 




Q 1 




1 n ] 






u, — 


1, - 


o 

^, - 


o 


4, - 


5, - 


ft 

D, - 


1, - 


8, - 


9, - 


10, - 


AT 

J-^P,hc,Ly,X 


2 


1 


1 


1 
























1 




















4 


4 


5 


8 


7 


3 


1 


















1 


5 


4 


2 
















36 


6 


25 


53 


56 


35 


17 


4 


1 














18 


49 


49 


34 


13 


4 












358 


8 


194 


454 


518 


404 


237 


100 


32 


6 


1 










158 


424 


488 


387 


220 


94 


26 


5 








3748 


10 


1590 


4036 


4953 


4324 


2947 


1565 


664 


208 


51 


7 


1 






1524 


3989 


4887 


4301 


2900 


1558 


641 


207 


44 


7 




40404 



TABLE VII. Factorization properties of characteristic polynomials CP{Tsq^Ly,di z) for cyclic strips of the 
square lattices. The notation {m,n) means that CP{Tsq^Ly,d, z) factorizes into polynomials in z of degree m 
and n with integer coefficients. The notation j'' indicates that there are k factors of degree j. 



Ly i d ^ 





1 


2 


3 


4 


5 


Np,sq,Ly,\ 


1 


(1) 


(1) 










2 


2 


(1) 


(1^) 


(1) 








4 


3 


(2) 


(1,3) 


(13) 


(1) 






10 


4 


(1,3) 


(4,5) 


(3,5) 


(1^,2) 


(1) 




26 


5 


(2.7) 


(8.13) 


(9.12) 


(5.8) 


(2.3) 


(1) 


70 
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TABLE VIII. Factorization properties of eharacteristie polynomials CP{Ttri,Ly.d, z) for cyelic strips of 
the triangular lattices. The notation (jn,n) means that CP{Ttri,Ly,d, z) factorizes into polynomials in z of 
degree m and n with integer coefiicients. The notation j'^ indicates that there are k factors of degree j. 



L,j I d ^ 





;i 


2 


3 


1 


•5 




1 


(1) 


(1) 










2 


2 


(1) 


(2) 


(1) 








4 


3 


(2) 


(1,3) 


(1,2) 


(1) 






10 


4 


(4) 


(9) 


(8) 


(4) 


(1) 




26 


5 


(2,7) 


(8,13) 


(9,12) 


(5,8) 


(2,3) 


(1) 


70 



TABLE IX. Factorization properties of characteristic polynomials CP{Thc,Ly,d^ z) for cyclic strips of the 
honeycomb lattices. The notation (m, n) means that the characteristic polynomial factorizes into polynomials 
in z of degree m and n with integer coefficients. The notation j'^ indicates that there are k factors of degree 
3- 



Ly i d ^ 





1 


2 


3 


4 


5 


^P,hc,Ly,\ 


2 


(1) 


(1^) 


(1) 








4 


3 


(3) 


(6) 


(1,3) 


(1) 






14 


4 


(1,5) 


(5,8) 


(4,7) 


(l^3) 


(1) 




36 


5 


(19) 


(43) 


(40) 


(22) 


(1,6) 


1 


132 
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TABLE X. Properties of the locus B for infinite-length cyclic/Mobius strips with width Ly of regular 
lattices A, including square (sq), triangular (tri), and honeycomb (he) lattices. The notation BCR denotes 
B crossings on the real q axis, the greatest of which is Qc- The notation SN indicates whether B has some 
support for values of q with Re{q) < 0, marked as y (yes) or n (no). The behavior of B in the neighborhood 
of the origin is indicated as convex or concave to the left, with the convention taken such that, for example, 
the circle |g — 1| = 1 is convex (to the left) at g = 0. The final column gives the reference. 



A 


Ly 


BCR 


SN 


B a.t q = 


Ref. 


sq 


1 


2, 


n 


convex 




sq 


2 


2, 


n 


convex 


[20] 


sq 


3 


2.34, 2, 


y 


concave 


[29] 


sq 


4 


2.49, 2, 


y 


concave 


[35] 


sq 


5 


2.58, 2, 


y 


concave 


[39] 


tri 


2 


3, 2, 


n 


convex 


[30] 


tri 


3 


3, 2, 


n 


convex 


[36] 


tri 


4 


3.23, 3, 2, 


y 


concave 


[36] 


tri 


5 


3.33, 3, 2, 


y 


concave 


here 


he 


2 


2, 


y 


concave 


[27] 


he 


3 


2, 


y 


concave 


[10] 


he 


4 


2.155, 2, 


y 


concave 


here 


he 


5 


2.26, 2.25, 2.20, 2, 


y 


concave 


here 
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XII. APPENDIX 3 



Here we list polynomials that occur in our explicit expressions for elements of transfer 
matrices. It is convenient to define 

One then sees that — 4g + 5 is ^4^3. Note that 

i^m+2,n. = (?-ir-'(g-2) (12.2) 

Fm+4,m = - ir-\q - 2)iq' -2q + 2) (12.3) 

Fm+6,m = {q- ir'-\q -2){q^-q + l)D^ 

- (? - ir-\q -2){q'-q + l){q' -3q + 3) (12.4) 
etc., so e.g., F4^2 — {q — ^){q — 2). One can also define 

Gm,n = - 2D,, (12.5) 

so that another frequently occurring entry, — 5g + 7 is seen to be ^4.3. We have not 
tried to carry out such constructions in an exhaustive manner and often just give shorthand 
notation to commonly occurring polynomials. Because of the various equivalences, there are 
also several ways of writing a given polynomial. 

We also list a set of polynomials that recur frequently and are given shorthand names. 



Pn, Tn, and Sn for quadratic, cubic, and quartic terms: 

P4 = g'-3g + 4, p6 = ?'-4g + 6, p8 = ?'-5? + 8 (12.6) 

Pio = - 6g + 10, pi3 = - 7g + 13, pu^q^ -7q + U (12.7) 

Pi5 = q^-7q + 15, pn = q^ - 8q + 17 (12.8) 

m = - 5q^ + Uq - 11, ns = - 6q^ + Uq - 13 (12.9) 

= - 7g2 + 19q - 19, rso = g^ - 7q^ + 19q - 20 (12.10) 
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= - V + 29g - 34, r46 = - Hg^ + 39g - 46 (12.11) 

= g^ - lOg^ + 36g - 47, = g^ - lOg^ + 37g - 49, rgg = g^ - llg^ + 43g - 58 

(12.12) 

s^^q"^ -5q^ + nq'^ -12q + 7, S21 = g'' - 7g^ + 21g^ - 32g + 21 (12.13) 
= g^ - 8g^ + 27g2 - 45g + 31, S47 = g^ - 9g^ + 34g2 - 63g + 47 (12.14) 

sg^ = g^ - 9g3 + 35g2 - 70g + 61, 8^2 = g^ - Hg^ + 48g2 - 99g + 82 (12.15) 
si48 = - 13g^ + 67g^ - 160g + 148, S173 = g^ - 13g^ + 68g^ - 171g + 173 (12.16) 
S3,5 = (?' - 2g + 3)(g^ - 4g + 5), t,,, = (g^ - 4g + 5)(g^ - 4g2 + 6g - 5) (12.17) 
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